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ABSTRACT 


This report is concerned with analytical investigation of si- 
nusoidally forced vibration of laminated, anisotropic plates including 
bending-stretching coupling, thickness-shear flexibility, all three 
types of inertia effects, and material damping. 

In the analysis the effects of thickness-shear deformation are 
considered by the use of a shear correction factor K, analogous to that 
used by Mindlin for homogeneous plates. Two entirely different approaches 
for calculating the thickness-shear factor for a laminate are presented. 
Numerical examples indicate that the value of K depends on the layer 
properties and the stacking sequence of the laminate. 

The general analysis is applicable to plates with any combination 
of natural boundary conditions at their edges. The analysis begins with 
the anisotropic stiffness and damping constitutive relations for a single 
layer and proceeds through the analogous relations for a laminate. Then 
the various types of energy and work terms are derived and the problem is 
formulated as an eigenvalue problem by application of the extended 
Rayleigh-Ritz method. 

The first five resonant frequencies of boron/epoxy plates with all 
edges free are calculated on the basis of two different resonance criteria: 
the peak-amplitude and modified Kenned y-Pancu techniques. The results show 
that the resonant frequencies obtained by the two techniques differ by only 
a very small amount, and are in good agreement with the results obtained 


iv 



both experimentally and analytically by Clary. Furthermore, the nodal 
patterns obtained agree satisfactorily with those of Clary. Finally, 
the damping values in this investigation are in good agreement with the 
experimental ones obtained by Clary. 
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SECTION I 


INTRODUCTION 


1.1 Introductory Remarks 

The continuing demand for increased structural efficiency in many 
advanced aerospace vehicles has resulted in development of laminated 
structures made of advanced fiber-reinforced composite materials, such as 
boron-epoxy, graphite-epoxy and boron-aluminum (reference 1). One of the 
most common structural elements for such vehicles is the rectangular plate 
or panel, for which numerous stable static, buckling, and free-vibrational 
analyses have been performed (reference 2). These advanced vehicles must 
maintain their structural integrity not only during high staticly applied 
mechanical and thermal loadings, but also in a variety of dynamic environ- 
ments. In order to predict the dynamic stresses to which a structure will 
be subjected, it is necessary to perform a dynamic response analysis in- 
cluding the effects of damping properties as well as the various stiffnesses. 
An example of this type of requirement is found in the Space Shuttle now 
under NASA preliminary development (reference 3) . 

A composite material is defined as two or more materials joined to- 
gether to form a nonhomogeneous material, which is used in constructing 
structures. Although such a material is nonhomogeneous on a micro scale, it 
behaves macroscopical ly as if it were a homogeneous, anisotropic material. 

An anisotropic material is one which exhibits different properties when 
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tested at different directional orientations within the body. For example, 
a single layer of composite material containing unidirectiona 1 ly oriented 
fibers (hereafter referred to as a unidirectional composite) has consider- 
ably greater stiffness and strength in the direction of the fibers than it 
does in a direction transverse to the fibers (see figure 1) . This aniso- 
tropic aspect is one of several which makes structural analyses for composite- 
material structures more complicated (reference 4) than those for ordinary 
isotropic materials, such as aluminum alloys or many unfilled plastics. 

A single layer of unidirectional composite has certain orthogonal 
axes of material symmetry and thus is said to be orthotropic . and behaves 
in considerably simpler fashion than a general anisotropic material. If 
the material is thin it is called plane orthotropic . However, for purposes of 
obtaining a more efficient design, it is often advantageous to place dif- 
ferent layers or plies at different orientations. For the case of a thin 
plate, this means that we must consider it as a plane anisotropic material. 
This treatment is intermediate between the complicated general anisotropic 
case and the much simpler plane orthotropic case. 

A laminate consists of two or more layers integrally joined together. 

It is said to he laminated symmetrica 1 ly if all of the layers above and be- 
low the midplane of the laminate have the same dimensions, properties, and 
|orientation| (if the layers are orthotropic); see figure 2. A laminate is 
said to be balanced if all of the layers oriented at +0 are balanced by an 
equal number of identical layers at an orientation of -0; see figure 3. 

If a laminate is not symmetrically laminated, coupling occurs between in- 
plane (either normal or shear) stress on one hand and either bending or 
twisting deformation on the other hand. In such a case as a laminate 
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consisting of a metallic substrate and a few layers of overlaying filamentary 
composite material, considerable bending-stretching coupling can occur and 
thus, it is necessary to consider not only the in-plane and bending stiff- 
nesses but also the coupling stiffnesses which couple together the other 
two effects (reference 4) . Obviously, this is considerably more complicated 
to analyze than a simple single layer. 

Unfortunately, it appears that no one has yet succeeded in devising 
a lamination scheme that is both symmetrical and balanced. However, in the 
'special, yet practical, case of multiple identical layers oriented alter- 
nately at +9 and -9, as the number of layers is increased, the bending- 
coupling effect diminishes (reference 5). For more than ten such layers, 
the bending-stretching coupling may be neglected for most engineering pur- 
poses and, thus, the laminate may be treated macroscopically as if it were 
a single-layer plate. 

It has long been recognized that, due to the relatively low shear 
stiffnesses of composite materials in planes normal to the laminating plane, 
the thickness-shear* deformations must be included in the analysis, even 
when the plate is relatively thin (references 6-8) , in order to achieve 
reasonably good predictions of laminate flexural behavior. Although 
several existing laminated plate theories include thickness-shear flexibility, 
all of them require an ad hoc assumption of the required shear correction 
factor. In this report, two different, relatively simple procedures are 

•jc 

Sometimes referred to as "transverse shear”, especially in the case 
of beams. The terminology used here follows that established by Yu for sand- 
wich plates (reference 9). Here the term, transverse, is reserved for the 
in-plane direction which is normal to the longitudinal direction. 
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presented to enable rational prediction of the appropriate composite shear 
correction factor (K) . The frequencies calculated using these values for 
K in shear-flexible laminated plate vibrational analysis give good agree- 
ment with the exact three-dimensional elasticity solution, which is com- 
putationally much more complicated. 

The use of high-damping polymeric materials in the form of a thin 
layer or tape has come into widespread use as a structural damper to reduce 
the vibrational response of aircraft panels, especially in high-noise regions 
such as in the vicinity of jet engines. When the polymeric material is 
added on either one side or both sides, it is known as an unconstrained 
damping layer; when the polymeric material is placed between two or more 
layers, it is called a constrained damping layer (or layers). Appropriate 
analyses and design procedures have been developed for including either type 
of damping layer, assuming all layers (metal and polymer) are isotropic 
(references 10-12). However, so far as known, no d eta i led solutions have 
been published on vibration of laminated composite-material plates which in- 
cluded material damping. Such an analysis is presented in this report and 
the results are compared with experimental data reported recently by Clary 
(reference 13) . 

1.2 A Brief Survey of Selected Vibrational 
Analyses of Laminated Plates 

Apparently the first vibrational analysis was carried out by Pister 
(reference 14) for a thin plate arbitrarily laminated of isotropic layers. 

In this case, the net effect of the bending-stretching coupling resulted in 
a reduced flexural stiffness. 
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Stavsky (reference 15) formulated a coupled bending-stretching 
dynamic theory for thin plates laminated of composite-material layers, 
but he did not present any numerical results. Apparently the first pub- 
lished results of the vibrational analysis of such plates is due to Ashton 
and Waddoups (reference 16) , who used the Rayleigh-Ritz method to analyze 
rectangular plates. These results compared reasonably well with experi- 
mental results for the completely free and cantilever cases. A similar 
analytical and experimental study, but involving simply supported and free 
boundary conditions, was carried out by Hikami (reference 17). Additional 
analysis of the free-edge case was carried out by Ashton (reference 18) . 

For simply supported plates, Whitney and Leissa (references 19,20) 
presented closed-form solutions for the natural frequencies in the case 
of cross-ply and angle-ply lamination schemes. As would be expected, their 
results showed a strong effect of bending-stretching coupling in lowering 
the natural frequencies. 

The case of clamped boundary conditions is more complicated analyti- 
cally, but more representative of practical aerospace structures. Rayleigh- 
Ritz and experimental investigations of such structures were carried out 
independently by Ashton and Anderson (reference 21) and by Bert and May- 
berry (reference 22) . 

The first more accurate vibrational analysis of laminated plates in- 
cluding thickness-shear flexibility was made by Ambartsumyan; see reference 
23. He assumed an arbitrary distribution of thickness-shear stresses through 
the thickness. However, in carrying out his actual calculations, Ambart- 
sumyan assumed a simple parabolic distribution. It can be shown by a 
simple mechanics-of -materials analysis (Jourawski shear theory; see Section 
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Ill) that the thickness-shear stress distribution must be discontinuous 
from layer to layer; thus, the simple parabolic distribution does not hold. 
Ambartsumyan did not give any numerical results for vibration of laminated 
plates; however, Whitney (reference 24) did so, using Ambartsumyan's basic 
theory. 

Using another approach, Yang et al (reference 25) entended the Mindlin 
homogeneous, isotropic, dynamic plate analysis (reference 26) to the laminated 
anisotropic case. They assumed a thickness-shear angle which is independent 
of the thickness coordinate (z) and then integrated the stress equations of 
motion to obtain the governing partial differential equations. After 
integration, they introduced a thickness-shear coefficient in an ad hoc 
fashion to correlate the predicted frequencies with known results. 

Also Yang et al introduced the coupling inertial effect (present 
only in the case of plates laminated unsymmetrical ly with respect to mass 
distribution), as well as the familiar translational and rotatory inertia 
effects present in homogeneous (or mass-symmetr ica lly laminated) plates. 

The inclusion of these higher-order inertial effects is consistent with the 
inclusion of thickness-shear flexibility. 

The analysis presented in this report is an improvement over both 
the Ambartsumyan and Yang et al analyses, in that it presents simple rational 
means for calculating the shear coefficient, rather than assuming it a priori . 

An alternative to considering shear deformation per se is to make a 
microlaminar analysis, such as considered by Biot (reference 27) and 
Bolotin (reference 28). However, this approach has not been used very ex- 
tensively so far. 

Another approximate approach is to use the Voigt and Reuss models to 
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determine the properties of an equivalent homogeneous, anisotropic, shear- 
flexible plate. This method was originated independently by Postma (refer- 
ence 29) , White and Angona (reference 30) , and Rytov (reference 31) to 
investigate wave propagation in a continuum consisting of alternating layers 
of stiff and flexible isotropic materials. This concept was applied re- 
cently to plates with experimental verification for the special case of a 
beam, by Achenbach and Zerbe (reference 32). A new analysis, extending the 
Postma approach to the more general case of a laminate consisting of alter- 
nating orthotropic layers, is presented in Appendix F. 

All of the analyses mentioned above are approximate formulations in 
this sense: interlaminar compatibility is impossible unless both of the 

in-plane Poisson's ratios are identical in all of the layers. The reason 
for this is that, in a plate, one can have discontinuities in the strain 
component in a given direction in the plane, due to the Poisson contraction 
caused by a stress resultant or a stress couple acting in the perpendicular 
direction in the plane. This deficiency can be removed by using the ap- 
proach introduced recently by Hsu and Wang (reference 33) and Wang (refer- 
ence 34) for laminated shells. Another technique is to apply the nonhomo- 
geneous three-dimensional elasticity approach, such as used recently by 
Srinivas et al (references 35-37) for a special case of simply supported 
edges. Still another method is to use finite elements in the thickness 
direction, as introduced recently by Tso et al (reference 38). Unfortunately, 
all of the more accurate analyses mentioned in this paragraph are quite 
complicated computationally and thus are not amendable to engineering 
analyses of practical structural elements. Furthermore, the work of refer- 
ence 37 showed that for structurally reasonable values of the following 
modal parameter, the Mindlin-type theory (such as extended here) is 
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sufficiently accurate for determination of natural frequencies, but not for 
determination of the associated stress distribution. 

There are numerous analyses in the literature on vibration of damped 
plates with isotropic layers; also there are a few vibrational analyses 
of single-layer, anisotropic plates and a few quasi-static (creep) analyses, 
of laminated, anisotropic plates. However, vibrational analyses of laminated, 
anisotropic plates are quite limited. Dong (reference 39) indicated the 
solution for the dynamic response of a simply-supported rectangular plate 
arbitrarily laminated of orthotropic, viscoelastic plies modeled as a 
standard linear solid. 
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SECTION II 


FORMULATION OF THE THEORY OF LAMINATED, SHEAR-FLEXIBLE PLATES 

In this section is presented a theoretical analysis of sinusoidally 
forced vibration of laminated, anisotropic plates including bending- 
stretching coupling, thickness-shear, all three types of inertia effects, 
and material damping. First the assumptions, on which the analysis is 
based, are stated explicitly. The general analysis is applicable to plates 
with any combination of natural boundary conditions at their edges. The 
analysis begins with the anisotropic stress-strain relations for a single 
layer and proceeds through the stiffness and damping constitutive relations 
formulation of the various types of energy and work terms, and culminates 
in the formulation of an eigenvalue problem by application of the extended 
Rayleigh-Ritz method. 

2.1 Hypotheses 

The following assumptions are made in the analysis presented here: 

HI. All displacements are assumed to be sufficiently small so that 
the linear strain-displacement relations are sufficiently accurate. 

H2. The layers which make up the plate are linearly elastic and 
may be isotropic or orthotropic with any arbitrary orientation in the 
plane of the plate. 

H3. The layers are sufficiently thin that thickness-normal-stress 
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effects may be neglected, i.e. the layers are assumed to have finite 
stiffnesses which resist membrane, bending, and thickness shearing 
stresses. 

H4. The layers are assumed to be bonded together perfectly. 

H5. The plate is assumed to have all components of translational, 
translational-rotatory coupling, and rotatory inertia. 

H6. All material damping effects are assumed to be small. They 
are incorporated by using stiffnesses which are complex rather than real 
(see Appendix B) ; the complex-stiffness array is assumed to be symmetric 
(see Appendix 1). All other thermal effects are neglected. 

H7. All initial-stress effects are neglected. 

H8. All interactions with a surrounding fluid can be either 
neglected or considered to be included in the values used for the material 
damping coefficients. 

H9. The excitation consists of a uniformly distributed normal 
pressure loading, sinusoidal with respect to time. 

2.2 Kinematics 

The displacement field is assumed to be (hypothesis H3) : 

u (x , y , z , t) = u q (x , y , t) + z^ (x ,y , t) 
v (x , y , z , t) * v q (x ,y,t) + z\|i y (x,y,t) (1) 

w (x , y , z , t) = w q (x , y , t) 

where u,v,w are the displacement components in the x,y,z directions (see 
figure 4) ; u 0 > v 0 » w 0 are displacement components at the middle surface 
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of the multi-layer plate; i|i and ill are the weighted middle-surface rotations 

x y 

in the respective xz and yz planes; and t is time. 

Within the framework of hypothesis HI, the total engineering strains 
are given by the following in-plane strain-displacement relations: 

(no sum) 

e ij = 6 i j + ZH ij (ij = xx »yy* x y) (2) 


where 


o 

e = u 
xx o,x 


0 

e ■ v 

yy °.y 



v + u 
o,x o, y 


(3) 


K = ill 
XX x,x 


X = ill 

yy y.y 


x’ = 111 + 111 

xy y ,x Y x , y 


(4) 


Also the following thickness-shear strain-displacement relations 

hold: 


G 


XZ 


W + ill 

o,x 


X 


'yz 


w + w 

o,y y 


(3) 


where a subscript comma denotes partial differentiation with respect to 

the variable following the comma, i.e. u == ^u /gx. 

o ,x o 


2.3 Stiffness Constitutive Relations 


In line with hypothesis H2 , the following stress-strain relations 
are assumed to hold for each individual layer of which the plate is com- 
prised : 
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I, 


f o (k) 1 

XX 


"" Q (k> 
H 11 

Q (k) 

W 12 

0 

0 

Q< k >“ 

W 16 


f (k) 

6 xx 

a (k > 

yy 


Q (k) 

^12 

Q (k) 

^22 

0 

0 

Q (k) 

V 26 


e (k) 

yy 

o< k > 

yz 

► = 

0 

0 


Q (k) 

W 45 

0 

< 

,< k > ► 
yz 

a (k) 

xz 


0 

0 

C 

Q (k> 

W 55 

0 


s (k) 

xz 

a (k > 

x y > 



Q< k > 

V 26 

0 

0 



(k) 

„ e xy J 


(k) (k) (k) (k) 

where a , a are the in-plane normal stresses; 0 .0 are the 

xx yy xz yz 

(k) 

thickness-shear stresses; a is the in-plane shear stress; the Q, . 

xy r ij 

are the reduced stiffness coefficients which are applicable to a thin 
layer (see Appendix A); superscript (k) denotes a typical layer k, where 
k = 1,2 , 3..., n ; and n is the total number of layers of which the plate 


( 6 ) 


is comprised. 

It is noted that it is more consistent mathematically to write 
the stress-strain relations in either matrix form as 


M 


(k) }= [Q^ (i,j=l,2,4,5,6). 


(7) 


or in tensor notation as 


_(k)_ (k) . 

°ij " Q ijkf e k£ (i,j-x,y,z). 


( 8 ) 


However, here a mixed notation is used for engineering convenience: 
double lettered subscripts (xx,yy,yz,xz,xy) following classical 
elasticity theory notation for stresses and strains, and double numbered 
subscripts (i, j=l ,2, 4, 5, 6) to reduce the number of subscripts used on 
the stiffness coefficients. 
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As is conventional in plate and shell theory, it is convenient to 

introduce generalized forces which are applicable to the whole laminate, 

rather than only a specific distance from the plate middle surface. 

(k) 

The stress components, , are functions of x,y,z; therefore, integrat- 

ing them through the thickness yields the following generalized forces, 
which are functions of x,y only: 


„h/2 


{n ,n ,n ,q ,q } = f {cr (k) ,o (k) ,o (k) ,o (k) ,o (k) ) ds 

( xx’ yy’ xy* x yJ J ' xx yy xy ’ xz yz J 


- f J (9) 

J l xx yy xy xz yz f 


k=l k-1 


»h/2 


{« ,M ,H ] . f {o (k) , c <k >, 0 <k >} 2 ^ 

I XX yy xyj J xx *y 

_ $ r Zk r (k) < k ) (k)i 

L J l. °xx ,a yy ,Q xy 1 2 d ‘ 


( 10 ) 


k=l ‘k-1 


where h = total plate thickness, the quantity (z^-z^ j) is the thickness 
of an individual layer (see figure 5), - membrane stress resultants 

(force per unit length) , =? thickness-shear stress resultants (force 

per unit length) , and are the stress couples (moment per unit length) . 

Inserting the stress-strain relations, given by equations (6) in- 
to equations (9) and (10), one obtains the following constitutive re- 
lations for the composite: 
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where the respective stretching, bending-stretching coupling, and 
bending stiffnesses of the plate are defined as follows: 


r^ 2 (k> 

Q ii 
-h/2 1J 


and the thickness-shear stiffnesses of the plate are as follows: 


r h/2 (k) 

; K, . Q. . dz 

u u J. h/2 u 


U,i-4,5> 


where = composite shear coefficient, introduced to account for the 
thickness-shear strain variation through the thickness. Section III 
presents several approximate analyses to predict K„ analytically. 

The A.j, > an d D.j submatrices appearing in equation (11) are 
present in the classical Kirchhof f-type theory of laminated plates (see 
refs. 4,5,15-22). In the case of an orthotropic plate all of the terms 
with subscripts 16 and 26, called the cross-elasticity or shear-coupling 


terms, vanish. Furthermore, when the layers are all isotropic materials, 
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only two of each set of the terms with subscripts 11,12,22, and 66 are 
independent . 

If the plate is either homogeneous (i.e. single layer) or laminated 
symmetrically about the plate middle surface, the bending-stretching cou- 
pling stiffness submatrix vanishes . Then the only remaining terms are 
the andD^j submatrices, which are present in classical, homogeneous 
thin-plate theory. 

The thickness-shear stiffness coefficients F. . account for the pre- 

1 1 

sence of thickness-shear strains, in a manner which can be considered to 
be a generalization of that used by Mindlin (ref. 26) for homogeneous, 
isotropic plates. To reduce the present theory to the classical Kirchhoff- 
type plate theory, the thickness-shear strains would be omitted; thus, 
from equations (5), one would obtain: 


Ilf * - w 

X 


o,x 


w 

°.y 


(14) 


Furthermore, one would delete the thickness-shear strains, e and e , 

xz yz 

from equations (11) and the thickness-shear stress resultants, Q and Q , 

x y 

would be computed from equilibrium considerations only. 


2.4 Strain Energy 


The differential of the strain-energy density (strain energy per 
unit volume) is given by: 


dU«. o <k) d e (k) + o (k) 
v xx xx yy 


de 


(k) 

yy 


+ a 


(k) 

xy 


de 


(k) 

xy 


15 



( 15 ) 


+a< k) d, (k) + o (k) de (k) 
xz xz yz e y2 

Integration of equation (15) yields: 


u< k) -U/2) |> (k) e (k) + o< k) e <« + o (k) e (k) 

yy yy xy xy 


xx xx 


(k) (k) L Jk) (k) 

J yz e yz 


+ a'"' e '"' + a e 

— xz xz J 


(16) 


The strain energy per unit of plate area (U ) is the integral of 

A 

00 

U over the total thickness of the plate: 
v 


n Z k 


U - f U^z-Y f u <k) dz 

A -h/2 v V , v 

k=l k-1 


(17) 


Substituting equations (16), (6), and (2) into equation (17), per- 
forming the integration, and using equations (12) and (13), one obtains 
the following expression for the strain energy per unit area: 


U. = (1/2) [A , (e° ) 2 + 2A,„ e° e° + A„(e° ) 2 + 2A., e° e° 

A 11 xx 12 xx yy 22 yy 16 xx xy 


+ 2A 0 , e e° + A,, (e ) + 2B. . g x 

26 yy xy 66 xy 11 xx xx 


+ 2 B 


1, - 

i ° ■ 

( e i 

k + e 

0 

x ) 

+ 2 B„ 

g° X 

12 

XX 

yy 


yy 

XX 

22 

yy yy 


(e° 

X 

+ 

o 

e 

X 

) + 2B 

( e° x + 

16 

XX 

xy 


xy 

XX 

26 

yy xy 

\ , r 

o 

£ 

X 

+ 

D, , 

CM 

+ 2 D„ 

xx + D, 

66 

xy 

xy 


11 

XX 

12 

xx yy 

I, f 

K 

X 

+ 

2 D 

K 

x + 

D X 2 +’] 

16 

XX 

xy 


/b 

yy xy 

66 yy 

» 

45 

e yz 

e xz 

+ 

*55 

2 

e xz 

] 



(18) 
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2.5 Damping Coefficients and Dissipative Energy 


For a material with Kimba 11-Lovell structural damping (see Appendix 
B) , the stress-strain rate relation for sinusoidal motion at a circular 
frequency U) can be expressed as follows: 


6 = (b/uu) e 


(19) 


where b is a material constant and e is the strain rate, 
The strain is given by 


e = lei e 


iwt 


( 20 ) 


where i = ,/~-l . 
Thus , 


e = ioue 


( 21 ) 


Substituting equation (21) into equation (19), one obtains: 


A * 

8 = Qe 


( 22 ) 


where 


Q = ib 


(23) 


Generalizing equation (23) to the entire array of stresses and 

strains in layer "k", in contracted notation analogous to equation (7), we 

* (k) 

obtain the following expression, where is symmetric (hypothesis H6) : 


K k) ( -tOM k) } 


(24) 
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Equation (24) is the same as equation (7), except that here the 
presence of the hat symbols ( *) denotes damping quantities rather 
than elastic quantities. 

The energy dissipated per unit of plate area, due to damping, 
is denoted by the symbol D^. The expression for it in terms of the 
midplane strains (e°j) and the curvatures (x.^) is the same as 
equation (18), except for the presence of the hat symbols over all 
of the A^, D ij‘ an< * quantities, where 


*tj * 1S 1J >A 1J • B 1J ' l8 iJ >B iJ ’ 


' l 8 M D lj • 


F s ig^F 
ij S i j ij 


(25) 


where g^ are loss tangents (see Appendix B) . 


2.6 Kinetic Energy 


The differential kinetic energy of an elemental volume (dx by dy 
by dz) is given by: 

(p (k) /2) (u 2 + v 2 + w 2 ) dx dy dz 
(k) 

where p is the density at the point (x,y,z) and u,v,w are the velocity 
components in the x,y,z directions respectively. 

Now the kinetic energy per unit of plate area is obtained from the 
differential quantity given above by dividing by dx dy and integrating 
over the entire thickness of the plate, as follows: 


18 



(26) 


(p (k) /2X<i 2 + v 2 + w 2 ) ds 


-h/2 


Substituting the kinematic relations, equations (1), into equation 
(26) and performing the integration, assuming that the density is uni- 
form through the thickness of each individual layer, one arrives at the 
following expression: 

T = (m /2)(u 2 + v 2 + vft + m. (u + v ’5/ ) 

A o o o o 1 o x o T y 


+ (m 2 /2) (^ 2 + ft 


(27) 


where m^m^, and m^ are respectively the mass per unit of plate area, 
first moment of mass per unit area (coupling inertia), and second 
moment of mass per unit area (mass moment of inertia) , defined as 
follows : 


{ m o ,m r m 2} 


t I 


.00 


i 1 , Z , z 2 ] 


dz 


(28) 


k=l 


k-l 


Only the coefficient m Q appears in classical Kirchhoff theory for 

the dynamics of thin, homogeneous plates. The Mindlin dynamic plate 

theory of homogeneous plates (ref. 26) contains both m^ and , while the 

dynamic theory originated by Yang et al (ref. 25) for laminated plates 

contains m , m , and m . 

o I 4 


2.7 Work Done by External Forces 
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The only external force considered here is a uniformly distributed 
normal pressure which has a sinusoidal wave form in time. Thus, it can 
be written in complex variable notation as follows: 

q(x,y,t) = q e^ (29) 

where q ~ amplitude of the normal pressure. 

The following expression for the total work done may be derived 
easily from the principle of virtual work: 

- =fr q (x.y.t) w q (x,y) dx dy (30) 

o o 


2.8 Application of the Extended Rayleigh-Ritz Method 

Since the present problem is one of steady-state harmonic ex- 
citation denoted in complex-variable exponential forms, the displacements 
and rotations are proportional to e* llt also. It is noted that the phase 
angle between the response and the excitation is taken care of by the 
imaginary components of these quantities. Thus, the time dependence 
cancels out in all of the energy terms which appear in Lagrange's equation. 
Therefore, it is necessary to consider only the amplitudes of the re- 
spective energy and work terms. Then, as shown in detail in Appendix 
C, the amplitude of the Lagrangian energy difference can be expressed 
as follows: 


L = (T + W) - ( U + D) 


(31) 
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where D,T,U,W are the amplitudes of the dissipative, kinetic, and strain 
energies and of the work done by the external forces, respectively. 

The strain and damping energy terms can be combined by introducing 
complex stif fnessses , as discussed in Appendix B. Making this sub- 
stitution and integrating over the plate area, we obtain the following: 
a b 

U+D = \ [ I {a., u 2 + 2A 10 u v + A 00 v^ + 2(A-./-u +A- v ) (v +u ) 

*’o o ^ 11 ° ,x 12 °>* o,y 22 o,y lo o,x 26 o>y o,x o,y 

+A,,(v +u ) 2 + 2B..u if + 2B 10 (u If +v iji )+2B v if 

66 o,x o,y 11 o,x T x,x 12 o,x T y,y o,y r x,x 22 o,y Y y,y 

+2B, [u (ip +if ) + (v +u ) i(i 1 + 2B 0 [v (if +i)i ) 

16 o,x Y y,x x,y o,x o,y x,x 26 o,y y,x x,y 

+ (v +u )lj ] + 2B (v +u )(i|i +if ) + 6 if 

o,x o , y y,y 66 o,x o,y Y y,x x,y ll T x,x 

+2D If If +6 „ if + 2D., If G +'lf ) +2D_ if (ft +'if ) 

12 Y x ,x Y y , y 22 Y y,y 16 x,x y,x x,y 26 Y y,y T y,x x,y 

+<jf ) + F. . (w +if ) + 2F.,(w, +ili ) (w +if ) 

66 v Y y , x Y x ,y 44 ojr y y' 45 ’y y o,x Y x 


+f 55 ( V + ^x )2 } dx dy 


(32) 


where the superscript (") denotes that the quantity is an amplitude, and 


V V 5 u ‘ <A i j +A i 3 ) • <W- <W ! M-i.2.6 


(33) 


f ij ‘ F tj +F lj ’ ‘’J" 4 ' 5 


The sum of the amplitudes of the kinetic energy and the work done 
by external forces can be expressed as follows: 
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( 34 ) 


T+W 






+ w ) + 2m, (u lit 
o 1 ox 



+ m 2 + i^)] dx dy + 


r j* 


q dx dy 


To apply the extended Raylelgh-Ritz method, the assumed functions 
for the amplitudes of the displacements and rotations are given the 
following general form: 


M N 


w 


- ) Y U $ (cr) * O) 

Li Li mn urn un 

m=l n=l 
M N 

= ) } V $ (a)* (p) 

L. L mn vm vn 

m=l n=l 

= / / W $(<*)♦ (p) 

' L mn wm wn 

m=l n=l 


M U 

/ / Y *. (cr)4> (P) 

/ : L ymn \Jjym fyn 

m=l n=l 
M N 

Y ) » *. (<y)*. (p) 

L L xmn (Jixm i|ixn 

m=l n=l 


(35) 


The modal coefficients U ,V ,W ,y are the undetermined 

mn mn mn ymn xmn 

parameters, the $'s are assumed modal functions, a ~ x/a, p y/b. 

Substituting equations (32), (34), and (35) into equation (31) 
for the Lagrangian energy-difference amplitude (L) , one obtains the 

result presented in Appendix Dl. 

Hamilton's principle can be stated mathematically as follows: 
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(36) 
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L dt 


= 0 


To achieve as close of an approximation as possible to equation 
(36) , the Lagrangian function L is minimized by setting its partial 
derivatives with respect to the modal coefficients respectively equal 
to zero, i.e. 


dL/*U k/ = 0 ; SL/dV^ = 0 ; dL/dW^ = 0 ; SL/SM' xk ^ = 0 


(37) 


^L/JVt' , = 0 ; k * 1,2,..., M ; f = 1,2,..., N 

yk-t 


Equations (37) represent a set of M x N nonhomogeneous , linear algebraic 
equations . 


2.9 Reduction to Matrix Form 


For computational convenience, it is desirable to express the set 
of equations (37) in matrix form as follows: 

^ tO - “■ [M] {5 J - { £ } < 3 « 

where [S] complex stiffness matrix, [M] ~ mass matrix, and ^? mn | and 
|f| are column matrices representing the generalized -displacements and 
generalized forces, respectively, i.e. 
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( 39 ) 


f \ 


r „ 

u 


0 

mn 



V 


0 

mn 



W 

III 

9 

mn 

i j 


V 


0 

' ymn 



Y 


0 

xmn > 




Equation (38) can be written in explicit complex form as follows: 


[s (R) + 


iS 


(I) 




^( R ) + 

• mn mn J 


- jj 


[m < r> ] { 5 < r > + IC 


'mn 


(I) 

mn 


} “ (* (W } 


(40) 


where superscripts (R) and (I) denote the real and imaginary parts of 
the complex quantities appearing in equation (38). 

Equation (40) can be solved for the response matrix, partitioned 
into real and imaginary parts, as follows: 


[«<*>] 

s< R >-A (R >' -S (I > 


t _ 


S (I > ! S (R> - A (R) 

( mn J 

t 


1 



J 


i 

f 



(41) 


The modal coefficients are placed in the inverse matrix in equation 
(41). For the specific numerical case treated M=N=2 , i.e. m,n=l,2. Thus, 
the complete matrix is of order 40 x 40 with each submatrix being of 
order 2x2. Thus, the problem has been reduced to a standard algebraic 
eigenvalue problem which can be solved by an available computer suboutline 
for the IBM 360 Series 60 digital computer. Complete computer program 
documentation is presented in Appendix J. 
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SECTION III 


DERIVATION OF THICKNESS -SHEAR FACTORS FOR LAMINATES 


In this section are presented two entirely different approaches 
for calculating the thickness-shear factors for a laminate. The first 
approach is to extend the Jourawski static theory of shear-flexible 
beams (reference 40) , as presented in elementary mechanics-of-materials 
textbooks, to a laminated beam. The other approach is to extend, to 
the laminate case, Mindlin's method (ref. 26) of matching the pure 
thickness-shear-mode frequencies predicted by two-dimensional dynamic 
elasticity and by Timoshenko one-dimensional, shear-flexible beam theory 
(references 41, 42). 

In the case of isotropic plates, there is only one independent 
thickness-shear factor (K) , since K^=0 and K^^*K^=K. The static 
approach mentioned above yields a value of K = 0.833 (reference 43) , 
while Mindlin's dynamic approach results in a value of 0.822 for K. 

Srinivas et al. (ref. 35) showed that use of either of these values for 
K in Mindlin's plate theory results in a very close approximation to 
the lower natural frequencies computed by exact, dynamic, three-dimensional 
theory of elasticity for rather thick, homogeneous isotropic rectangular 
plates simply supported on all edges. 

In the case of orthotropic plates, there are two independent thick- 
ness-shear factors (K, . and K cc ) , since K. c =0 and K__^K. .. In the case 

44 55 45 55 44 

of a plane anisotropic plate, such as one consisting of a unidirectional 
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composite material with its major material-symmetry axis oriented at an 
acute angle (9) with its edges, there are three independent, non-zero 
However, for this case, Appendix A presents transformations from 
which the three F^, proportional to K„ and defined in equation (13), 
can be calculated from the following data: F and F for the orthotropic 

case (9=0°) and the value of 9 in the anisotropic case. Thus, the general 
problem of determining is reduced to one of calculating and K,.,. 
for the orthotropic case. 

To calculate K,.,. for an orthotropic laminate by either of the two 

approaches mentioned above, one considers a beam oriented in the x 

(k)* (k) 

direction and laminated of layers having properties E^ and C,.,. . To 


calculate K /(/) , x -• y 


(k) (k) r (k)^ r (k) 

> E 22 * and C 55 C 44 * 


To provide a check for the results, the fundamental eigenvalues 
for free vibration calculated by laminated, shear-flexible plate theory 
(Section II), using these values of K^, are compared with the exact 
laminated elasticity theory values given by Srinivas et al. (ref. 35). 


3.1 Static Approach 


To calculate the horizontal shear force per unit width, AF, acting 
on a cross section, figures 6 and 7 are used. For consistency plate notation 
is used, even though the theory is only one dimensional. 

The bending stress in typical layer "k" , at a distance z from the 
midplane of the laminate, is calculated as follows: 


(k) c (k) 

o = E ' x z 

XX 1 1 XX 


(42) 


* (k) 

For thin plates as treated in this report, - 


«(k) . 

Q n , where 


0 (k)_ E (k), (k) (kK 

Q 1 1 * E 11 ' U v 12 V 21 ; • 
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where x^ == longitudinal bending curvature and h longitudinal Young's 

modulus of layer "k". 

The longitudinal stress couple (bending moment per unit width) is 
defined as follows: 


M 

xx 


n Z, 

y r V k > 

/.< - XX 


z dz 


k=l z 


k-1 


(43) 


The longitudinal curvature can be expressed in terms of M x as 


follows : 


X = M /D.. 

XX XX 1 1 


(44) 


where =• longitudinal flexural rigidity. 

From equations (42) and (44) , the bending stresses acting on the 
left (x=x^) and right (x = x^+Ax) sides of the elements are: 


* 'C 2 /D n> M xx ; 




(43) 


Since both of these bending stress distributions act on identical 
cross-sectional areas, the horizontal shear force per unit width, which 
acts on the element, is given by: 


= J [ °xx )(x l^ x >-°xx )(x l )] dA 
A 

o 


where A =? area shown in figure 8 . 
o 


(46) 
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Combining equations (45) and (46), one obtains: 


AF = (AM xx /D ll )Y(z) 


( 47 ) 


where 


Y(z) = J E^z dA 


11 


The force per unit width, AF, must equilibrate the horizontal 
(k) 

shear stress, o , acting on the bottom face of the element shown in 
xz 

figure 7. Thus, we have: 


AF = Ax 

xz 


(48) 


Hooke's law in shear for layer "k" is 


a (k) . c (k) £ (k) 
xz 55 xz 


(49) 


( k) ( k) 

where and e are the modulus and strain corresponding to longi- 

55 xz 

tudinal thickness-shear action in layer "k". 

Substituting equation (49) into equation (48) , one obtains the 
following result: 


AF = C^ k) e (k) Ax 
55 xz 


(50) 


The following expression for the longitudinal shear strain at 
a distance z is obtained by equating the right-hand sides of equation 
(47) and (50) : 
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(51) 


e^(x,z) = [Y(z)/(C^ } D u )] Q x (x) 

where is the thickness-shear stress resultant (force per unit width), 
which is related to M^^by static equilibrium as follows: 

Q = AM /Ax (52) 

X XX 

The shear strain energy for a laminated, rectangular-cross-section 
beam, shown in figure 9, is: 


U . **J. C 5? ts « l 2dA < 53 > 

A 

The shear strain energy associated with an equivalent uniform, 
longitudinal thickness-shear strain in a laminated, rectangular-cross 
section beam of cross-sectional area A is: 


U' = 
s 



dA 


(54) 


where K,._ is the longitudinal thickness -shear factor and e' is weighted 
55 xz 

according to: 


xz 


clj 

k=l : 


,(k) 


. (k) 

a e 
xz xz 


d ^ [ l J d *T‘ (55 > 


k-1 


k=l 


'k-1 


Equating the right-hand sides of equations (53) and (54) yields an 
equation; then substituting equations (49) and (51) into that equation, 
one arrives at the following explicit expression for the longitudinal 
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thickness-shear factor 


n z, 


c £ J Y( Z ) dz y 


k=l z 


k- 1 


55 


n z, 


^ Z k 


( 56 ) 


[ Y J dz] ( y I ( c “)'‘w«)f dz} 


k=l k-l 


k=l k-l 


Equation (56) is a general expression for the static, longitudinal 

thickness-shear factor for an arbitrarily laminated beam. It is a 

relatively simple algebraic expression which depends upon only the 

longitudinal Young’s and thickness-shear moduli of the individual layers 
(k) (k) 

(Ejj and ' ) and the lamination geometry. Numerical results are 
presented in later examples. 


3.2 Dynamic Approach 

Due to the inherent mathematical form of the equations of dynamic 
elasticity theory for a multi-material (multi-layer) medium, it is not 
feasible to present a general solution for the thickness-shear factor 
obtained by using the dynamic approach. Even in the case of a specific 
class of laminate, such as the three-layer, symmetrically laminated one, 
it is not possible to obtain an explicit expression for the thickness- 
shear factor. In fact, the complexity of the dynamic case, which in- 
volves the dynamic elasticity-theory analysis of a shear-flexible 
laminated beam, approaches that of the shear-flexible laminated plate and 
thus negates much of the advantage of using the thickness-shear factor 
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approach to laminated, shear-flexible plate dynamics. However, in order 
to illustrate the dynamic approach and to provide a numerical comparison 
with the static approach and with the exact results of Srinivas et al 
(ref. 35), an analysis for the symmetrical three-layer case is presented 
in Appendix E. 
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SECTION IV 


MODAL FUNCTIONS USED FOR VARIOUS 
PLATE BOUNDARY CONDITIONS 


In applying the extended Rayleigh-Ritz method, the approximate 
modal shapes assumed must satisfy the kinematic boundary conditions. 
However, to improve convergence, it is desirable to satisfy the force- 
type boundary conditions also. For a rectangular plate, the three cases 
of boundary conditions most commonly encountered in practice are discussed 
in Sections 4. 1-4.3. 

It should be noted tha C , C , Z , Z are characteristic parameters 

m n m n r 

tabulated in reference 44. Also, five boundary conditions per edge need 
to be prescribed rather than four as in classical theory. 

4.1 Simply Supported on All Edges 


The boundary conditions considered here are simply-supported 
in the sense that 


N 
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= 0 

at x = 0 ,a 
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yy 

yy ° 

0 X 




edges 


(37) 


As was pointed out by Wang (reference 45) , it is impossible for a 
separable-form deflection function to satisfy the above boundary conditions 
for the case of a plate of anisotropic material (Q and Q„, ^ 0). However, 
the following modal functions permit eqs. (35) to satisfy eq. (57) exactly 
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in the orthotropic case, but only approximately in the anisotropic case: 
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For the simply supported case, we substitute equations (58) into 
equations (D-l through D-6) , The evaluations of the integral forms are 
presented in Appendix G. 


4.2 Fully Clamped on All Edges 

The boundary conditions for a fully clamped laminated plate are as 
follows : 
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The assumed modal functions selected to permit equations (35) to 
satisfy equations (59) are as follows: 


$ = $ = $ =4 = s in 2mrro 

um vm tyxm ijiym 


$ = $ =$, =$, = sin 2nrr8 

un vn ijixn ijiyn 


(60) 

$ = cosh Z a - cos Z o - C (sinh Z a - sin Z a) 

wm m m m m m 

$ = cosh Z p - cos Z 6 - C (sinh Z p - sin Z P) 

wn n n n n n 


33 



For the fully clamped case, we substitute equations (60) into 
equations (D1 through D6) . The evaluations of the integral forms are 
presented in Appendix G. 


4.3 Free on All Edges 

The boundary conditions appropriate for a laminated plate free on 
all edges can be expressed mathematically as follows: 


N = M =0 at x = 0 , a and y = 0 , b 

xy xy 
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(61) 
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(D1 through D6) . The evaluation of the resulting integral forms are 
presented in Appendix G. 
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SECTION V 


NUMERICAL RESULTS AND COMPARISON WITH 
RESULTS OF OTHER INVESTIGATORS 


In this section a comparison is made between the numerical results ob- 
tained in this investigation and those obtained by other investigators. 

5.1 Thickness-Shear Factors for Laminates 


Two example problems are analyzed and, where possible, results of 
Example 1 are compared with solutions obtained from other sources. The 
results of Example 2 are presented for use in design. 

Example 1: Three-layer, symmetric, isotropic laminate . - The coordinate 

system is shown in figure 10. All geometrical and material property ratios 
correspond to the data of ref. 35 given as follows: 


2 /z x = h (2) /[h (2) +2h (i) ] = 0.8 , = 15, 

v = 0.3 


XI)/. (2) 


<D/o<2) = ,c .(1)/J2) =1 . (1) = (2) _ 


C 55'/ C 55'“ 15 ’ P /P =1 . and v 


Assuming only the x-direction is considered in the static and dynamic 
approaches which were developed in Section III .then the following results 
are obtained: (I) static case - the result obtained from equation (56) 

can be expressed as K^=0.651, (II) dynamic case - the result obtained by 
equating eqs. (E-13) and (E-27) is K<,^=0.350. 
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When the three-layer, laminate thickness parameter is chosen to be 
z 2^ z l > t * ie analysis given in Section III can be reduced to a single- 
layer homogeneous, isotropic material. Thus, using the data E^/E^ = 
^55^55^ = P^/p^^ = l anc * v (1) = v (2) = 0.3, the static case gives 
K,.,- = 0.833 and dynamic case results in K,.,. = 0.822. These results are 
identical with the values of the static and the dynamic shear factors 
obtained for single-layer plates by previous investigators (refs. 43 and 
16) . 

Example 2: Multiple alternating layers of two materials . - The 

geometrical and material property ratios are given as 


h (2)/ h ( 1 ) „ = C^/C^ 2) = from 0.01 to 100 


The laminates vary from two layers to nine layers. Since the static 

case is a relatively simple algebraic expression which considers only 

longitudinal Young's and thickness-shear moduli of the individual layers 
(k) (k) 

(E^ and )and the laminate geometry, for design purposes, eq. (56) 

will be used. The results are shown in table 1 and figure 11. 


5.2 Plate Simply Supported on All Edges 

Here we consider only free vibration of plates made of composite 
material having material-symmetry axes coinciding with the geometric 
coordinates (plate edges) as shown in figure 12 . Neglecting energy 
dissipation, the time integral of the difference between the potential and 
kinetic energies attains a stationary value. The displacement and rotations 
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are proportional to e 1Ujt . The time dependence in the strain and kinetic 
energies cancel each other, i.e. only the amplitudes of the strain energy 
and kinetic energy need be considered. Thus, we use only the real part 
of homogeneous eqs. (D1-D6) with proper boundary conditions, see eq. 

(57) , to apply Example 3. 

Example 3: Three-layer, symmetric , isotropic plate . - The elastic 

coefficients are those of isotropic material: 


q<» . qOO . ,<«/{, . [v <»] 2 } >Q <« - v (k) E <k) /{l-[v (k> ] 2 } 
Q^ k) - q* k) * Q* k) - E (k) / [2[l+v (k) ] j , q‘ k >-0 (k-1,2) 

Geometrical and material property ratios correspond to the data 
of ref. 35, given as follows 

z 2 /z 1 = h (2) /[h (2) +2h (1) ]= 0.8, eJ^/eJ 2 ^ C^/C^-15, 


(1) . (2) _ (1) _ (2) 

p / o =1, v = v =0.3 


and the following dimensionless geometric parameter defined in ref. 35: 
2 I r ... _i, (k) i2 r ^(k) /. -)2 _ 


0.0002n 2 for k=l 


g 2 = | [mnh (k) /a] 2 + [nnh (k Vb] 2 ={ „ 

k ; l 0 . 0128 tt for k=2 


When the thickness of three-ply laminates is chosen to be z 2 /z^ = 1, 
it can be reduced to a single-layer material as a special case. Then the 
material properties are as follows: 
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i-.j = 1> 2, 4, 5, 6 




!<»/ E<f - c£/C<?> = p<‘>/ p (2) - . ; v<» - v (2 > - 0.3 


g l = °» s 2 = 12 * 


A FORTRAN IV program Is employed to compute the lowest eigenvalues 
of plates with simply supported edges. The two cases considered are: 

(1) single- layer homogeneous, isotropic, and (2) three-ply symmetrical 
construction. The results are shown as a solid line in figures 13 and 
14. The static and dynamic shear factors calculated in examples 1 and 
2 are also shown for comparison. 

The static Jourawski shear theory and dynamic Timoskenko type 
theory appear to give good results for K as evaluated by inserting in 
laminated, shear flexible plate theory and comparing the lowest eigen- 
values with those given in refs. 35 and 14. 

5.3 Plate Free on All Edges 

The problem considered here is forced vibration of a rectangular 
plate made of an arbitrary number of composite-material layers each 
having its major material-symmetry axis oriented at an arbitrary angle 
with geometrical coordinates as shown in figure A-2. 

Since the present problem is one of steady-state harmonic excitation 
including material damping effects, the complete set of eqs. (D1-D6) with 
proper boundary conditions, eqs. (61), must be used. 
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Example 4: Boron/epoxy plate with twenty-four parallel plies 


oriented at an arbitrary angle. - The input-data material properties are 
Young's and shear moduli (E^, ^22* ^44* ^55* ^66^’ ma j° r Poisson's ratio 
(v^) , bending and twisting* stiffnesses (D^, ^22* ^12’ ^66^’ anc * t * ieir 
corresponding loss tangents, as generated in ref. 46 from constituent- 
material experimental properties. It is noticed that all of the above input 
data are functions of frequency. Geometrical and material properties 
correspond to the data of ref. 13, as follows: 


length, a = 18.19 in ; width, b = 2.75 in. 

thickness, h = 0.034 in. ; density, p = 0.000194 lb-sec /in 
angle of orientation, 0 = 0°, 10°, 30°, 45°, 60°, 90° 

A FORTRAN IV program is employed to compute the eigenvectors of 
a plate with all edges free. To obtain the resonant frequencies of the 
first five modes, two different resonance criteria are applied. In the 
peak-amplitude method, the system is excited harmonically and the amplitude 
at a particular point (o = 0 , g = 1) is measured over a range of fre- 
quency. The amplitude of the response depends not only on the dynamic 
characteristics of the system, but also on the amplitude of the force 
applied to it. For a linear system the peak amplitude is taken as 
maximum displacement per unit amplitude of force. Throughout this re- 
port, this ratio is referred to as the response. A typical result of 
response with varying frequency is shown in figure 15 for an angle of 
orientation, 0 = 10°. The complete results are as tabulated in Table II. 

In the modified Kennedy-Pancu method, the modal shapes, phase 
relations between motions at various points, and coupling between the 


various degrees of freedom are assumed to be unaffected by the presence 
of a small amount of damping. Vectors are used for the analysis of the 
results of forced vibration by making a polar plot (called the Argand 
plane) of the displacement vector (amplitude and phase angle) . In the 
vicinity of resonance, the resulting curve would be an approximately 
circular arc. Each point on the circle locus corresponds to a value of 
frequency. 

In the modified Kennedy- Pancu method, the resonant frequency is de- 
fined to correspond to the point of the Argand-plot curve having maximum 
change in arc length (As) per fixed change in frequency (Aw) . The 
change in arc length As can be determined approximately in analytical 
fashion as the response amplitude (W/q) multiplied by the difference in 
the phase angles corresponding to the two distinct frequencies (separated 
by Aw) , i . e . 

As = (W/q) Acp 

The parameter As/Aw is plotted versus excitation frequency in figure 
16 for a specific angle of orientation, 9 = 10°. 

As can be seen in the figure, the curves exhibits sharp spikes and 
thus the Kennedy-Pancu resonant frequencies, i.e. the frequencies cor- 
responding to the maximum points, can readily be found. The results obtained 
in this fashion for all of the plates used by Clary are tabulated in Table II. 

A discussion of experimental methods used to determine material 
damping in composite materials is presented in Appendix H. Here the damping 
ratios are calculated by using the modified Kennedy-Pancu method as expressed 
in equation (B-69), Appendix B. It is noted that the damping ratio is the 


40 



ratio of material damping coefficient to critical material damping coef- 
ficient. The results associated with each mode are tabulated completely 
in Table II. 

The nodal patterns of the first five modes are calculated from the 
response at selected points on the surface of the plate. Twenty-five 
points in the x direction and five points in the y direction are used. 
The nodal patterns are shown in figures 17-22. The results appear to 
be in good agreement with the data of ref. 13. 
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VI. CONCLUSIONS 


The analyses in Sections II and IV were developed for the vibra- 
tional problem of composite-material plates with thickness-shear flexibility 
and material damping. The theory used in the analyses is that of Yang, 
Norris, and Stavsky (ref. 25), which can be considered to be the laminated, 
anisotropic version of Mindlin's dynamic plate theory. 

To account for the effects of thickness-shear deformation, a shear 
correction factor K was introduced. A variety of methods to arrive at 
an appropriate shear factor have been proposed for homogeneous plates; 
the most popular ones resulted in K = 0.833 for a static distribution 
(ref. 44) and 0.822 for the dynamic case (ref. 26). Srinivas et al. 

(ref. 35) showed that use of either of these values for K in Mindlin's 
plate theory gives a very close approximation to the lower natural fre- 
quencies computed by exact, dynamic, three-dimensional elasticity theory 
for rather thick, homogeneous, isotropic rectangular plates simply sup- 
ported on all edges. 

Unfortunately to date no one has proposed a means of rational 
calculation of the shear factor for a laminated plate. Srinivas et al. 
made an exact, dynamic three-dimensional elasticity analysis of simply 
supported laminated plates. However, their analysis is quite tedious 
computationally and thus they have presented numerical results for only a 
very limited number of cases. 
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In the present Investigation, both static (Jourauski shear theory) 
and dynamic (Mindlin type analysis for pure thickness-shear motion) 
approaches were used to derive a shear factor K for laminates. 

An assessment of the accuracy of these theories was made by com- 
paring the lowest eigenvalues calculated by the laminated, shear-flexible 
plate theory with the laminated, three-dimensional exact values given 
by Srivinvas et al. (ref. 35). 

It has been found that the shear factor K for a laminate is not 
the same value as that for a homogeneous, isotropic member. The numberical 
examples indicated that the value of K depends on the layer properties 
and the stacking sequence of the laminate. 

Using peak-amplitude and modified Kennedy-Pancu methods, the first 
five resonant frequencies were calculated for various laminated boron/ 
epoxy plates with all edges free. The resonant frequencies obtained by 
the two techniques differed by only a very small amount, and were in good 
agreement with the results obtained both experimentally and analytically 
by Clary (ref. 13). Furthermore, the nodal patterns have been found to be 
satisfactory or as good as the data given by Clary. Finally, the damping- 
ratio results were in good agreement with the experimental ones obtained 
by Clary. No comparison with analytical results for lmainated-plate 
damping could be made, since they have not been available previously. 
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APPENDIX A 


NOTATION AND TRANSFORMATION FOR ELASTIC COEFFICIENTS 

The generalized Hooke's law (Cauchy equations) for a general 
anisotropic material in terms of rectangular coordinates x,y,z are as 
follows: 
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where 0 . . and g.. are stress and strain components, and C c are the 
ij ij kl 

Cauchy stiffness coefficients for a three-dimensional body. It can 
be shown that the matrix of Cauchy coefficients is symmetric, provided 
that the material is conservative (elastic) , so that 


C tk = 


(A-2) 


Usually a single layer of fiber-reinforced composite material has 
the fibers oriented more or less parallel to the top and bottom surfaces 



of the layer, i.e. material-symmetry plane XY coincides with surface 
plane xy. Thus, it is orthotropic in the yz and xz planes, and the 


shear-normal or cross-elasticity coefficients with subscripts i4 and 
i5 (i=l,2,3,6) vanish: 
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The only exception to the above, i.e. the only case when the general 
form (A-l) , rather than (A-3) , must be used is the case of a shingle- 
laminated composite, as shown in figure A-l (reference 47). 

If, in addition to being parallel to the top and bottom faces of 
the layer, the fibers are also oriented in a direction parallel to one 
of the coordinate directions in the plane of the layer, the material- 
symmetry plates XY, YZ ,ZX all coincide with the coordinate planes xy, 
yz.zx and the material is said to be completely orthotropic . In this 
case the Cauchy relations are: 
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(A-4) 



where the asterisks denote the orthotropic case. 

In the case of a thin layer of fiber-reinforced composite 
material in which the thickness -normal stress (c^^can be neglected 
(see Hypothesis H3, Section 2.1), the third equation in set (A-4) can 
be solved for e^ z in terms of and e^. Then equations (A-4) can be 
simplified as follows: 



(A-5) 


where the reduced orthotropic stiffness coefficients (Q^) are related 

to the Cauchy three-dimensional orthotropic stiffness coefficients 
* 

(C^) as follows: 
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(ij = 44,55,66) 


Equation (A-5) can be written in the form of a single matrix 
equation as follows: 

l°IJ i = { e Ij} 

where 
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Qn Q 12 0 0 0 
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44 

0 0 2Q* 5 0 

0 0 0 Ke 


(A-9) 


where the presence of the factor 1/2 in the e^j matrix is used to make 
it a second-rank tensor and this necessitates the presence of the factor 
2 appearing in equation. (A-9) . 

In structural-panel applications of composite materials, usually there 
are design requirements for multiple orientations of fiber-reinforced com- 
posite-material layers. Therefore, it is essential that a set of trans- 
formation relations be used to calculi the stiffness coefficients for 
any desired orientation from that associated with the major material- 
symmetry direction (fiber direction) as shown in figure A-2. Such 
relations are developed in the ensuing paragraphs. 
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Arbitrary orthogonal axes in the plane of the laminate and making 


an angle 9 with the materia 1-symmetry axes (X,Y) are designated as 
the x,y axes. Then the components of stress and strain can be trans- 
formed from the x,y axes to the X,Y axes as follows: 
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C V { «ij } 


(A- 10) 
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and 


where and are as defined in equations (A-8) ; 

Kj} are similar except that the subscripts X,Y,Z are replaced by 
the subscripts x,y,z; and [T 1 is the transformation, defined as 


follows : 
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m - cos 9 , n sin 9. 


For the case when 9 t 0, we have 



In view of equations (A- 7 ,A-10 ,A-11 ,A-14) , 
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By substituting equations (A-9) and (A-12) and then performing 
the matrix operations indicated in equation (A-15) the elements of 
the [Q matrix are related to those of the [Q ki ] matrix as follows 
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APPENDIX B 


MODELS AND MEASURES OF MATERIAL DAMPING 
Bl. Mathematical Models for Material Damping 

Internal friction or damping in materials can be caused by a variety 
of combinations of fundamental physical mechanisms, depending upon the 
specific material. For metals, these mechanisms include thermoelasticity 
on both the micro and macro scales, grain boundary viscosity, point- 
defect relaxations, eddy-current effects, stress-induced ordering, inter- 
stitial inpurities, and electronic effects (reference 48) . According to 
Lazan (ref. 48), little is known about the physical micromechanisms 
operative in most nonmetallic materials. However, for one important class 
of these, namely polymers and elastomers, considerable phenomenological 
data have been obtained. Due to the long-range molecular order associated 
with their giant molecules, polymers exhibit rheological behavior inter- 
mediate between that of a crystalline solid and a simple liquid. Of 
particular importance is the marked dependence of both stiffness and damp- 
ing on frequency and temperature. 

The purpose of developing a mathematical model for the rheological 
behavior of a solid is to permit realistic results to be obtained from 
mathematical analyses of complicated structures under various conditions, 
such as sinusoidal, random, and transient loading. According to reference 
48, as early as 1784 Coulomb recognized that the mechanisms of damping 
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operative at low stresses may be different than those at high stresses. 

Even today, after more than 2500 publications on damping have appeared, 
major emphasis is placed on linear models of damping for several reasons: 
they have sufficient accuracy for the low-stress regime and linear analyses 
are computationally more economical than nonlinear ones. 

The simplest mathematical models of rheological systems are single- 
parameter models: (1) an idealized spring, which exhibits a restoring 
force linearly proportional to displacement and thus displays no damp- 
ing whatsoever, and (2) an idealized dashpot, which produces a force 
linearly proportional to velocity. Obviously, neither of these models 
are very appropriate to represent the behavior of most real materials. 

The next most complicated models of rheological systems are the 
two-parameter models: (1) the Maxwell model, which consists of a spring 

and dashpot in mechanical series, as shown schematically in figure B-l(a) ; 
and (2) the Kelvin-Voigt model, which is comprised of a spring in 
parallel with a dashpot, as shown in f igure B-l (b) . The Maxwell model 
is a fair approximation to the behavior of a viscoelastic liquid. How- 
ever, as a model for a viscoelastic solid, it has several very serious 
drawbacks (reference 48); there are no means to provide for internal 
stress and for afterworking. The Kelvin-Voigt model overccmes these 
deficiencies and is a first approximation to the behavior of a viscoelastic 
solid. However, it has the following disadvantages (reference 48): there 
is no elastic response during application or release of loading, the 
creep rate approaches zero for long durations of loading, and there is 
no permanent set irrespective of the loading history. 
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It should be mentioned that the Kelvin-Voigt model is the simplest 
one which permits representation as a complex quantity when subjected to 
sinusoidal motion. To show this, one can write the following differential 
equation governing the motion of a single-degree-of-freedom system con- 
sisting of a sinusoidally excited mass attached to a Kelvin-Voigt element: 


mu+cu+ku=Fe 


iO)t 


(B-l) 


where m 2 mass, c 2 viscous damping coefficient, k 2 spring rate, F s 

exciting force amplitude, uj b circular frequency of exciting force, 

u 2 displacement, i = V - 1 , t s time, and a dot denotes a derivative 

with respect to time. Assuming that sufficient time has passed for the 

* 

transients to die out, we can write the steady state solution of equation 
(B-l) as follows: 


~ i(wt-cp) ~ 2 . . , -.-1 i&t 

u = u e = F [(k-muu )+ iiucj e (B- 2 ) 

where cp 2 phase angle between the response (u) and the exciting force, 
and 11 a displacement amplitude. 

It is noted that the terms containing Kelvin-Voigt coefficients k and 
iuc can be grouped into a single Kelvin-Voigt complex stiffness (k) as 
follows : 


*It can be shown that they do die out for a mechanical system, for which 
c, m, and k are all positive quantities. 
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k = k + i'juc 


(B-3) 


Alternatively, the Kelvin-Voigt element coefficient could have 
been grouped into a single complex damping coefficient (c) as follows: 


c = c - (i k/u,) 


(B-4) 


The force (F^) in a dashpot is given by 


= c u = ioucu e 


i (out -c$ 


(B-5) 


The energy dissipated per cycle (U^) is a dashpot undergoing 
sinusoidal motion can be calculated as follows: 


,.2n/jLi 


U d = F^du = J F^ u dt = ncu/. u 


(B-6) 


In 1927, Kimball and Lovell (reference 49) observed that many 

engineering materials exhibited energy losses which were contradictory 

to equation (B-6) . Specifically, they found that was proportional 
~2 

to u but independent of jj, i.e. 


U d = C i U 


( B— 7) 


~2 

Later work by Wegel and Walther (reference 50)also showed that 1^ ^ u 

but C' was a weak function of ju. Still later, for stresses below the 
d 

~2 

fatigue limit of the material, Lazan (references 48 ,51) showed that U^u 
and was practically independent of x. 
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Equating the energy dissipated per cycle in the Kimball-Lovell 
material, equation (B-7), to that in the dashpot, equation (B- 6 ), one obtains 

C ^ = ttcuj (B- 8 ) 

To utilize the Kimball-Lovell observation in practical structural 
dynamic analyses, it was desirable to have an expression for the as- 
sociated damping force, i.e. an equation analogous to equation (B-5) , which 
governs a dashpot. According to Bishop (reference 52) this was first done 
by Collar. Combining equations (B-5 , B- 8 ), one obtains the following 
"frequency-dependent damping" relation: 

F d = (b/uu) u (B-9) 


where b is a constant given by 

b = C'/tt (B-10) 

a 

The kind of damping represented by equat ion (B-9) has been called fre- 
quency-dependent damping, since the usual dashpot coefficient c is now 
replaced by the quantity b/uj. 

For the damper represented by Eq.(B-9),it is convenient to replace 
k in the undamped system by the following Kimball-Lovell complex stiffness 
(reference 53 ) : 

k' = k + ib (B-ll) 

where it is now noted that both k and b are assumed to be independent of 
frequency. This representation has been used extensively in aircraft 
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structural dynamic and flutter analyses (references 54,55). 

An alternative to equa tion (B-ll), is to use the concept of a complex 
damping coefficient, originated by Myklestad (reference 56), in which the 
spring constant is replaced by 

k* ' = C e im (B-12) 

where and m are constants. 

Although Kimball and Lovell's original work was based on data 
obtained during steady-state sinusoidal motion, the models represented 
by equations (B-9, B-ll, B-12) have been applied to free vibrations as 
well (refs. 52,53,57). To alleviate the difficulty of the interpreting of 
uj in equation (B-9)for the case of free vibration or for multiple-fre- 
quency forced vibration, Reid (reference 58) suggested the following form: 

F d = b | u/u | u (B-13) 

where the bars denote that the absolute value of the quantity between 
them is to be used. 

In summary, four different versions of material damping have been 
discussed. In terms of the differential equation for free vibration, 
they may be written as follows: 


mu' 

+ (b /■'!,) u + ku - 0 

(B- 14a) 

mu 

+ (k + ib) u = 0 

( B- 14b) 

mu' 

+ Cj e im u - 0 

(B-14c) 

mu* 

+ bju/u|u + ku“0 

(B- 14d) 
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The models represented by equations (B-14) have been criticized 
by many investigators for various reasons: 

1. In equation (B-14a) , the interpretation of a, is dubious (refer- 
ences 58,60). Milne (reference 59) suggested that it is associated with the 
imaginary part of the pair of complex roots, but went on to state that 
there is no clear justification for this interpretation. 

2. Equations (B-14b) and (B-14c) are equations with complex coef- 
ficients and the meaning of this is not clear (references 57 , 59) . Further- 
more, although they have complex exponential solutions, neither the real 
nor imaginary parts alone are solutions (references 59,60). 

3. Equation (B-14d) is a nonlinear differential equation, due to 
the behavior of |u/u] (reference 59). 

4. None of the models represented by equations (B-14) can be sim- 
ulated, even on a conceptual basis, on an analog computer (references 

58, 60). 

5. The model represented by equation (B-14a) does not meet the 
Wiener-Paley condition (references 61,62) of causality for physically 
realizable systems, as was shown in references 60, 63-65.. 

6. Equation (B-14b) fails to give the proper relations of a 
damping force that is proportional to displacement and in phase with 
velocity. This can be alleviated by use of a more unwieldy expression 
proposed in references 66,67. 

Incidentally reference 68 claimed that none of the equation 
(B-14) models result in energy dissipation that is independent of 
u; (the reason for abandoning the Kelvin-Voigt model). However, an error 
in this reasoning, which makes it invalid, was pointed out in reference 
69. 
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It can be shown that the sinusoidally forced vibration version of 


free vibration equations (B-14a,b,c) are exactly equivalent. Probably 
the most popular form for writing this is as follows: 

mu' + (k + ib)u=F e l (B-15) 


The major criticism of the use of equation (B-15) is that it results 
in a conventionally defined magnification factor (see Section B2) , correspond- 
ing to zero frequency, which is not equal to unity for g / 0 (reference 52) . 
This difficulty can be eliminated by redefining the magnification factor as 
follows: 

MF' = u/u' (B-16) 

st 

where 



T/k' 


and k' is the stiffness coefficient associated with the total in-plane force 
amplitude, i.e. 


k' = (k 2 + b 2 )* 


Thus , 


MF' = (k 2 +bV/[(k-mu; 2 ) 2 + b 2 ]* = (1 + g 2 )V[(l-w 2 /uj 2 ) 2 +g 2 ]^ 

n 


The major advantages of using the equation (B-15) model for steady - 
state vibrational and flutter analyses are as follows: 

1. The analysis is considerably simplified in comparison with that 
for viscous damping, as pointed out in refer® ce 70. 

2. The analysis is linear. 

3. The complex -modulus concept has long been in common use, especially 
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in connection with measurement of the damping properties of elastomers, 
polymers, etc. (reference 69) . 

4. The complex-stiffness approach has been used very successfully 
in flutter analysis for approximately three decades refs. 54, 55, 69, 71). 

In view of these considerations, the complex-stiffness approach is 
used in the main portion of the present report. However, for completeness 
and comparison, a number of other more complicated approaches are also dis- 
cussed in this section. 

To overcome some of the deficiencies of the Kelvin-Voigt and Max- 
well models, they were combined in various ways to obtain the three-para - 
meter models shown in figure B-l(c) and (d) . It can be shown (reference 
72) that by properly selecting the values of the coefficients, the two 
models can be used interchangably . In other words, the model repre- 
sentation is not unique. In fact, both representations have been called 
the "standard linear solid", "standard linear material," "simple an- 
elastic model", or "standard model of a viscoelastic body" (references 
73-76). Unfortunately, few materials have been characterized by the use 
of the standard linear solid. Perhaps the most extensive characterization 
has been carried out for flexural vibration of 2024-T4 aluminum alloy: 
by Granick and Stern (reference 77,78 ) in air and by Gustafson et al 
(reference 79) in vacuo. 

Obviously one can go on from a three- parameter model to a four- 
parameter one, etc. Continuing this process, one obtains two well- 
known models (reference 74): the Kelvin chain, shown in figure B-2(a), 

and the generalized Maxwell model , shown in figure B-2(b). The be- 
havior of such systems can be written in either differential or integral 
form, as discussed in detail in reference 74. However, the complexity of 
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such a model has limited its use both experimentally (in characterizing 
engineering materials) and analytically (in performing dynamic analyses). 
Thus, for engineering applications, there is a continuing search for 
models which are simple, yet represent the behavior of real materials 
with sufficient accuracy. 

A ^relatively simple model, which is somewhat similar to the generalized 
Maxwell one, is due to Biot (reference 80 ) ; see figure B-2(c) . The number of 
simple Voigt elements in this model is allowed to increase without bound, 
so that the damping force is represented by the following expression: 



(B-17) 


where t is a dummy variable; e,g^= parameters, end Ei(u) is the exponential 
integral 


Ei (u) 5 J (e‘^/?)d5 (B-18) 

CO 

Caughey (ref. 60) made a detailed study of Biot's model for both 
free and forced vibration. He showed that for a/e > 10, the energy dis- 
sipated in Biot's model is within 4% of being independent of frequency (a) 
and thus essentially depends only upon the square of the displacement 
amplitude. Milne (ref. 59) studied the Biot model by means of the con- 
volution integral. 

Neubert (reference 81) introduced a variant of the Biot model, in 
which a different distribution function for the local stiffness-damping 
ratio is used. Milne (ref. 59) also studied this model via the convolution 
integral. Milne introduced a new synthesized model of his own. However, 


59 


the main conclusion of his investigation was that, for practical engineer- 
ing purposes, the choice of a rheological model is not critical, provided 
that the range of constant damping is kept the same and the damping is 
sma 1 1 . 

The linear hereditary theory of material damping was originated in 
1876 by Boltzmann (reference 82). In this theory, the energy loss is 
attributed to the elastic delay by which the deformation lags behind 
the applied force. It is called a hereditary theory because the in- 
stantaneous deformation depends upon all of the stresses applied to the 
body previously as well as upon the stress at that instant. The damping 
force is given by 

t 

F d “ j §(t,r) u(r) dT (B-19) 

-09 

where t == actual time, t = an instant of time (dummy variable), u = 
displacement, and $ = hereditary kernel or memory function. Often the 
hereditary kernel depends upon the difference (t-T) only; then equation 
B-19) becomes: 


F 


d 


j f(t - t) u (t) 

V 


dT 


(B-20) 


It is most advantageous to determine the hereditary kernel from 
experimental data. It has been found to be a monotonical ly decreasing 
function which can be represented mathematically as follows: 


n 

i (t) = £ A.e' ait 
i=l 


(B-21) 


Such a system is said to have a heredity of degree n and A^ and q. 
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are called heredity constants. Hobbs (ref. 75) has shown that a first- 
degree hereditary model is equivalent to a simple anelastic model 
(standard linear solid). 

Volterra (reference 83) has presented a method of determining the 
hereditary kernel from material response to either suddenly applied 
loading or to sinusoidal excitation. He also presented plots of mag- 
nification factor and phase shift verus frequency for a wide range of 
parameters of a single-degree-of-treedom system with first-degree 
hereditary damping. For certain combinations of the parameters, the 
magnification and phase characteristics are quite different from that 
of a viscously damped system; yet for certain other combinations, the 
characteristics are quite similar. 

Numerous nonlinear mathematical models have been proposed to per- 
mit better correspondence between theory and experiment. However, all 
of them, by their very nonlinear nature, are more complicated to apply 
in structural dynamic analyses of practical engineering systems. Perhaps 
the simplest one is damping energy per unit volume and per cycle pro- 
portional to stress to a power (ref. 48): 

U dv = C d o m (B-22) 

where and m are material constants and <x is the cyclic stress 
amplitude. It is noted that m=2 corresponds to the Kimba 11-Lovell re- 
lation (linear theory) discussed previously. It has been found that m 
varies from 1.8 to 6.0, generally being near 2.0 at low stress amplitudes. 

The hysteresis loop associated with the complex-modulus model, 
equation (B- 11), is elliptic in shape; see figure B-3(a). However, 
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the hysteresis loops determined experimentally have sharp corners at both 
ends and are nearly linear in the low-stress, low-strain region, as shown 
in figure B-3(b) . According to ref. 48 (page 97), in 1938 Davidenkov 
proposed a nonlinear mathematical model which results in a similar shape 
of hysteresis loop. Further work on this class of model has been carried 
out by Pisarenko (reference 84) . 

Rosenblueth and Herrera (reference 35) proposed a nonlinear model 
which eliminates the causality difficulties of the complex -modulus model. 
Chang and Bieber (reference 86 ) introduced a nonlinear hysteresis model 
in which there is no hysteresis damping unless the displacement amplitude 
exceeds a certain threshold value. 


B2. Measures of Material Damping 

In this section, various definitions of damping are discussed in 
the context of a complex-modulus material. 

Energy Dissipation Under Steady-State Sinusoidal Vibration . - Energy 
dissipated per cycle (U^) in the form of internal frictional heating is 
one measure of damping. However, this quantity depends upon the size, shape, 
and dynamic stress distribution (which in turn, depends upon the particular 
mode of vibration). In view of the above difficulties, the specific damping 
energy (U^^) is usually considered to be a more basic property of the 
material . rather than the structure. The is defined as the damping energy 
per cycle and per unit volume, assuming a uniform dynamic stress distribution 
throughout the volume considered. Thus, the total damping energy is 


U 


r 


IT dV 
d v 


(B-23) 
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where V : volume. The usual units of U, are in-lb/cycle and of U. , 

d dv 

3 

in-lb/in -cycle. 

Resonant Magnification Factor Under Steady-State Sinusoidal Ex - 
citation. - The resonant magnification factor (RMF) is defined as the 
dimensionless ratio of the response at resonance to the excitation, 
where both the response and the excitation must be specified in the 
same units (see figure B-4) . These quantities may be in units of dis- 
placement, velocity, accelerat ion , or strain. Unfortunately , the re- 
sonant magnification factor is dependent upon the structural system con- 
figuration as well as upon the damping property of the material, so that 
it is considered to be a system characteristic, rather than a basic material 
property. 

Incidentally RMF is not applicable to nonlinear systems, since then 
it depends upon the excitation level. (In a linear system, RMF is independ- 
ent of the level of excitation.) 

Bandwidth of Half-Power Points Under Steady-State Sinusoidal Excitation . 
The separation (u^-uUj) between the frequencies associated with the half- 
power points increases with an increase in damping (see figure B-5) and thus 
can be used as a measure of damping. A more meaningful definition is to 
normalize by the associated resonant frequency (uu n ) , i.e. use the dimension - 
less bandwidth given by: 


(X 2 ' J 'l )/J 'n 

This measure is the basis for determination of damping by the original 
Kennedy-Pancu method (ref erence 87 ) . 

The quality factor Q is defined as the reciprocal of the dimensionless 
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bandwidth . 


Derivative of Phase Angle with Respect to Frequency, at Resonance . 
This measure is the basis for determination of damping by an improvement 
(reference 88) of the Kennedy-Pancu method. 

Loss Tangent Under Steady-State Sinusoidal Excitation . - Applying 
the complex-stiffness (see Section Bl) to the material property represent 
ing stiffness, namely the Young's modulus (E) , we obtain 

E - E R (1 + ig) (B-24) 

ic 

Thus, the loss tangent is defined as follows: 

g ? E l /E R (B-25) 

J 'k'k'k 

where E = loss modulus and E is called the storage modulus . 

Referring to figure B-6, we obtain the following relation: 

g = tan Y = E I /E R (B-26) 


or 


v = arctan E*/E R 


(B-2 7) 


*Also known as the "loss coefficient", "loss factor", or "damping 
factor". 

**Also sometimes called the "dissipation modulus . " 


***Also sometimes called " elastic modulus" or "real modulus." 



The quantity y is called the loss angle and equation (B-26) explains 


the origin of the term loss tangent for g. 

Cyclic Decay of Free Vibrations . - In a linear system, free vibrations 
decay exponentially (see figure B-7) . The larger the damping, the faster 
is the decay. Thus, the logarithmic decrement 6 is defined as follows: 

6 = In (a i /a i+1 ) (B-28) 

where In ? natural logarithm. 

For a power-law material, eq. (B-22) , to have 6 independent of 
amplitude, m must be 2 (see Section B3 for proof). When m=2, the follow- 
ing means of calculating 6 is more practical, especially when damping 
is small: 

6 = (1/n) >.n (a./a ) (B-29) 

l i+n 

where n is any arbitrary integer. 

Temporal Decay of Free Vibrations . - Another measure associated with 
the decay of free vibrations is the temporal decay constant v^, defined as 
follows (ref erence 48 , foldout opposite p. 35): 

v t - (^-t^ 1 >.n [u(t L ) /u(t 2 )] (B-30) 

where t^ and t ^ are two different values of time, arbitrary except that 
t 2 > t^ and (t 2 - t^) must be an integer number of periods of damped 
vibration. The quantities u(tp and u(t 2 > are the respective displacements 
at times t^ and t 2 * It is seen that the usual units for v t are sec \ 
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An alternate way of specifying the temporal decay is the deca y 
rate y^, which has units of decibals per second (db/sec) and is defined 
as follows: 

Y t = 20 1 log C u ( t 1 ) /u( t 2 ) ] (B-31) 

or 

-1 V t^ t 2 _t P 

Y t = 20(t 2 -t 1 ) log e L = (20 log e) v » 8.68 v 

where log = logarithm to the base 10. 

Spatial Attenuation of a Plane Wave in a Slender Bar . - If 
one propagates a plane, harmonic wave in a slender bar made of a homo- 
geneous, linear material, the wave exhibits exponential decay with axial 
position (x) , analogous to the temporal decay of free vibration (which 
may be consiered as a standing wave). Analogous to the logarithmic de- 
crement (see above), we have the logarithmic attenuation 6 , defined as 
follows in terms of amplitudes a(x^) and a(x 2 ) at stations x^ and x 2 : 

6 g = tn [a(x 1 ) /a(x 2 ) ] (B-32) 

Analogous to the temporal decay constant, we define the spatial 
attenuation constant v as follows: 

s 

V 5 (x -x ) 1 In [a(x,)/a(x,)] (B-33) 

S 4 1 1 4 

The unit of v is in l 
s 

Finally, analogous to the decay rate, we define the spatia 1 

attenuation rate v as follows: 

s 

i 
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(B-34) 


Y g - 20 ( x 2 _x p 1 log [a(x 1 > /a(x 2 > ] ~ 8.68 v g . 

The unit of y is db per unit length. 

It should be cautioned that the determination of damping by wave 
attenuation requires that the cross section of the bar be compact (pre- 
ferably round or square) and that the largest cross-sectional dimension 
be very small compared to the wave length of the traveling wave. Further- 
more, the two stations along the bar, at which the wave measurements are 
made, must be: (1) sufficiently far from the point of impact that initial 

transients have died out and (2) sufficiently far from each other so that 
a measurable change in amplitude occurs. Kolsky discussed these points 
in detail in reference 89 . 

B3 . Inter-Relationships Among Various Measures of Damping 
for Homogeneous Materials 

Here we consider a single degree-of-f reedom system consisting of a 
mass supported on a massless spring made of a homogeneous Kimba 11-Lovell 
or complex-modulus material. Then the governing differential equation 
for simple harmonic excitation is equation (B-15), which has the following 
steady-state solution: 


where 


u = u e i(j ' t-c P) 

(B-35) 

~ F [(k-mu, 2 ) 2 +b 2 ]’ l/2 

(B-36) 

2 

cp s arctan [b/(k-m&> )] 

(B-37) 


The magnification factor (MF) is defined as follows: 
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where u is the so-called static displacement, defined as follows: 

S L 

u = F/k (B-39) 

S L 

Combining equations (B-36, B-38, and B-39), one obtains the 
following relation: 


MF = {[1 - (uj/uu ) 2 ] 2 + (b/muu 2 ) 2 }" l/2 (B-40) 

l n n J 

where uu n is the natural frequency (i.e. the resonant frequency of the 
undamped system) , given by: 

1 /2 

d, n « (k/m) ' (B-41) 

The critical material damping coefficient (b c > is the minimum value 
of b for which the free vibratory motion is non-oscillatory. It 
is obtained from the following relationship: 

(b /2mo.' ) 2 = k/m = uu 2 
c n n 

or 

b = 2ma. 2 (B-42) 

c n 

The damping ratio £ is defined as follows: 

Q - b/b c (B-43) 

Combining equations (B-40, B-42, and B-43), one arrives at the 
following general, dimensionless relationship: 



MF = {[l-(tu/uJ n ) 2 ] 2 + 4£ 2 }" 1/2 


(B-44) 


It is seen easily that MF is a function of only two parameters: the fre- 
quency ratio and the damping ratio (0 . It should be mentioned 

that equation (B-44) is somewhat different than the known relationship 
for a single-degree-of-freedom Kelvin-Voigt viscously damped system, 
even though the damping ratio £ is defined analogously ( s c/c , 

where c = 2muj ) . 
c n 

The resonant magnification factor (RMF) is defined as the value 
of MF at resonance, which is taken here to occur when cp = 90°. In view 
of equation (B-37) , this implies that the resonant frequency is uo^. 

Thus, at resonance, equation (B-44) reduces to the following simple 
expression: 


RMF = (20 


(B-45) 


Thus, it is seen that the same simple relationship between the damping 
ratio and the resonant magnification factor holds as does for a Kelvin- 
Voigt system. Although the resonant amplitude can be measured accurately, 
there are problems in measuring the static deflection in complicated 
structures (due to signa I-to-noise ratio limitations of instrumentation). 
For this reason, in structural dynamics, RMF is seldom used as a measure 
of damping. 

The damping energy, i.e. the energy dissipated per cycle, can be 
computed as follows: 

, f.2n/vi; 

= <j> F^ du = j F^iidt (B-46) 


where F^ is the damping force, given by 
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F, = ( b / di) u 
d 


(B-47) 


Substituting equations (B-35) and (B-47) into equation (B-46) and 
integrating, we obtain the following result: 

IT = rrb u 2 (B-48) 

d 

To convert equation (B-48) to the form of equation (B-22) , with 
m=2 , it Is necessary to divide by the volume and to relate u to a. 

Now the force amplitude is 

F = A o = ku (B-49) 

where A - cross-sectional area. 

Assuming the spring is a massless and uniform bar, we have 

k = AE/L (B-50) 

where L 5 effective length of spring and E - Young's modulus of the bar. 
Combining equations (B-48, B-49, B-50), one obtains the following result; 

U. = U./V = (nbL/A) (o/E) 2 (B-51) 

d s d 

Comparing the form of equations (B-22) and (B-51) , we see that 
for a Kimba 11-Lovell material, m = 2 and 

C = rrbL/AE 2 (B-52) 

d 

It should be noted that Lazan (reference 48 , chap. II) has shown 
that these results are independent of the stress distribution and the 
specimen geometry, provided that m = 2. 
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In view of equations (B-42, B-43, and B-44), equation (B-52) can 


be expressed in the following more useful form: 


C d = 2tt C/E 


(B-53) 


According to equation (B-48) , the damping energy (and thus, 
excitation power) is proportional to the square of the displacement 
amplitude u q (and thus, the MF) . Then the power is one half of what it 
is at resonance when MF = RMF//2. In view of equation (B-45) , the half- 
power magnification factor (HPMF) is related to the damping ratio as 
follows : 


HPMF = RMF//2 = (2/2 f) 


(B-54) 


To determine the two sideband frequencies (u^ = and cor- 
responding to half power, we combine equations (B-44) and (B-54) to 

2 

obtain the following quadratic expression in u; : 


( 2/2 C ) 2 = [1 - (u>/u , n ) 2 ] 2 + 4 Q 2 


which has the following positive roots: 


^ ( 1±2 0 k (B-55) 

1,/ n 

It is interesting to note that equation (B-55) is identical to the analogous 

2 

expression for the lightly damped Kelvin-Voigt system, i.e. ( c/c ) « 1. 

For small damping (f « h) , 



Thus, we get the following relation between the dimensionless bandwidth 
and the damping ratio Q: 


(u^-ujj) / uj^ ss 2£ 


(B-56) 


In the original Kenned y-Pancu method (reference 87) , this relationship is 
used to determine the damping ratio from the geometry of an experimentally 
determined Argand plot (a polar-coordinate plot of response amplitude 
versus phase angle) . 

In terms of the quality factor Q, equation (B-56) may be rewritten 
as follows: 


Q ~ (20 


(B-57) 


It is interesting to note that in order for the half-power frequencies 
to be real, the damping ratio £ must be no greater than /2/2 « 0.707. How- 
ever, this is not a severe limitation for actual structural materials. 

The potential energy (strain energy) stored in our single-degree-of - 
freedom mathematical model is given by: 



du 


(B-58) 


where F g is the spring force, given by 

F = ku 
s 


(B-59) 


Thus , 

u 

U = k J u du = ku^/2 (B-60) 

o 

Combining equations (B-48) and (B-60) , we obtain 


U /U = 2-nb/k 
d 


(B-61) 





Using equations (B-41, B-42, B-43), we can rewrite equation (B-61) 
in the following more useful form: 


or 


U./U = 4tt£ 

a 


(B-62) 


C = (l/4n)(U d /U) 


(B-63) 


The loss tangent is defined as follows for the system considered 

here : 


S = b/k (B-64) 

Inserting equations (B-41, B-42, B-43) into equation (B-64) we 
obtain the following useful relationship: 



(B-65) 


The following useful relationship, which is sometimes taken as a 
fundamental definition of damping (reference 66), is obtained here by 
combining equations (B-63) and (B-65) : 


g = ( l/2rr) (U /U) 
a 


(B-66) 


Taking the derivative of phase angle cp with respect to frequency 
di, one obtains the following result from equation (B-37) : 


dcp/djj = 2bnW[ (k-tna;^) ^+b^ ] 


(B-67) 


At resonance (subscript R) , defined by = /k/m, equation (B-67) 


becomes : 


(dcp/du;) = 2mu, /b 
K. n 


(B-68) 


Combining equations (B-42, B-43, and B-68), we obtain the following 
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very useful result : 


C = [w n (dqp/daj) R ] 


(B-69) 


It Is Interesting to note that although the MF vs. ui relation, 
equation (B-44) , is different than the corresponding relation for a Kelvin- 
Voigt viscously damped system, equation (B-69) is identical to its cor- 
responding relationship for a Kelvin-Voigt system. Equation (B-69) is used 
to determine the damping ratio from an experimentally determined Argand 
diagram in the improved Kennedy- Pancu method (reference 88). 

The transient solution of the system represented by equation (B-14a) 
can be written as follows: 


u' = "O'e , ’ a n t sin CAl-f 2 oj t - «p) 

n 


(B-70) 


From its definition in equation (B-28) , the logarithmic decrement 6 
can be determined as follows: 


"O t 
J n 1 


6 = (n 


ruu t 

’ n 


(t ~ T) = ln e = (B-71) 

* n 1 


n 


where T is the period of damped oscillation, given by 


2 - 1/2 

T = (2n/ JJn )(l-C ) 


(B-72) 


Thus , we now have 


6 = 2tt C (1-CV 1/2 


(B-73) 


For small damping (^ «1) , equation (B-73) can be replaced by the 
following approximate expression: 


6 ~ 2n Q 


(B-74) 
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Combining equations (B-65) and (B-74) , one obtains the following: 


6 = TTg 


(B-75) 


In view of equation (B-63) and since the specific energies are 
given by: 


we have 


u d S ■ V v - u 8 * u/v - 


r = U /4nU 
ds s 


However, the specific strain energy is 


(B-76) 


(B — 7 7) 


U = a / 2E 
s 


(B-78) 


A power- law-damping material is defined by equation (B-22) . 
Combining this relation with equation (B-78), one obtains: 


Q = (E C d /8n)o' 


m-2 


(B-79) 


Thus, it is clear that Q is independent of stress amplitude only for a 
KL (Kimball-Lovell) material (m=2) . Then 


C « E C d /8n 


(B-80) 


or, in view of equation (B-65) , 


g = E C d /4rr 


(B-81) 


Furthermore, since the logarithmic decrement of a KL material with 
small damping is approximately proportional to the damping ratio as shown 
by equation (B-74), the following approximate relation holds: 


6 « E CV4 

d 


(B-82) 
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APPENDIX C 


DERIVATION OF ENERGY DIFFERENCE 


It is assumed that the generalized coordinate § and generalized 
force f have the following forms: 


S = 5 e 


f = q e 


iujt 


iajt 


(C-l) 


where | is a complex form representing u ,v ,w , it and It . 

r r °ooox T y 

The strain energy can be expressed as 


n n 




2iout 

U e 


(C-2) 


i=l j=l 


where Q are stiffness coefficients and U is the amplitude of the strain 
energy given by: 


n n 


i=l j=l 


(C-3) 


For the damping energy, a dissipation function suitable for material 
damping (reference 90) is developed. Hence, we introduce a dissipation 
energy function for material damping as follows: 


D = 


M 11 

*1 l 

i“l j=l 


(b ij /uj) 


l t Sj =(W2) 


n 

r 

L 

i“l 


n 

Y A 
L" ij 
j-1 


n _ «*• 2iu;t , . 2iuut 

Q, i e =(iU))D e 


(C-4) 
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where b. . are structural damping coefficients, Q, . = ib. . and D is the 
ij ij ij 

amplitude of dissipation energy given by: 


n n 


* ' * 1 l \i hh 

i=l j=l 


(C-5) 


The kinetic energy of the system is 


n n 


n n 


T = % Y y m. . . = - (oj 2 /2> Y y m. .?.?• e 2lUJt =-T e 

L L ij *i J L L i;ri*j 

i=l j=l i=l j=l 

(C-6) 


2iujt 


where m^ are inertia coefficients and T is the amplitude of kinetic 
energy given by 


n n 

T = (J/2) Y Y m i . ?. %. 

i=l j=l 


(C-7) 


The work done by the uniformly distributed normal force is 


W = f y F =i = q e iUJt y f. - w e 
i=l i=l 


2iuut 


(C-8) 


where q is the amplitude of the normal pressure and W is the amplitude of 
the work done 


n 

w - q L S t 
i=l 


(C-9) 


The Lagrangian equation takes the following form (reference 91) 


d ,*1* • SD 
dt ( ar + a£ 



(C-10) 


Each term of equation (C-10) can be expressed in terms of its 
amplitude as follows: 



3T iuut 
— e 

35 


bd 


BD iuut 
— e 

5? 


(C-ll) 



35 


35 


iuut 

e 


With equations (C-ll) 
(C-10) as follows: 


, f 


aw _ aw iuut 
35 35 e 


introduced, one can rewrite equation 


||- (T + W - U - D) e iUJt = 0 (C-12) 

or 

(T+W-U-D) =0 (C-13) 

35 

Equation (C-13) is analogous to one presented by Volterra (reference 83 ) . 

Thus, the amplitude of the Lagrangian energy difference can be 
expressed as follows 

L = (T + W) - (U + D) (C- 14) 


when used in conjunction with Hamilton's principle, expressed in equation 
(36). 



APPENDIX D 


COMPLETE ENERGY EXPRESSIONS 


Dl. The Energy Difference 


M N 11 


V r f (tA /a 2 )(U 4’ 4 ) 2 +(2A 19 /ab)(U 4' 4 )(V 4 4' ) 

L L J J L 11 mn um un 12 mn um un mn vm vn 

m=l n=l 0 0 

+2A 1 .(a“ 1 U 4' 4 )(a _1 V 4' 4 +b~ 1 U 4 4' ) 

16 mn um un mn vm vn mn um un 


. ^ 


+2A_.(b -1 V 4 4' )(a _1 V 4' 4 +b" 1 U 4 4' )+A,,(a _1 V 4' 4 +b ~ lu mi A, m $ ' > 

26 mn vm vn mn vm vn mn um un 66 mn vm vn mn um un 

+2B 11 (a" 1 U 4' 4 )(a _1 Y 4' 4, )+2§ 1 -[(a" 1 U 4' 4 )(b' 1 'V 4 4' 

11 mn um un xmn ijfxm tyxn 12 mn um un ymn tyym tyyn 

+b _1 V 4 4' )(a _1 Y 4', 4, )]+2B. 0 (b _1 V 4 4' )(b‘V 4, 4', ) 

mn vm vn xmn ijfxm \|/xn 22 mn vm vn ymn ijjym ,fyn 

+2B 1 ,L(a" 1 U 4' 4 ) (a _1 v 4.' 4* +b' 1 'f 4. 4.' )+a _1 V 4' 4 

16 mn um un ymn y ym \jlyn xmn ijixm \|)xn mn vm vn 

+b _1 U 4 4’ )(a _1 Y 4' 4, )]+2B 0 ,[b' 1 V 4 4’ )(a‘ 1 l' 4,’ 4, 

mn um un xmn Y xm fxn 26 mn vm vn ymn \jj ym yn 

+b" 1 Y 4, 4! )+(a _1 V 4' 4 +b _1 U 4 4' Kb' 1 ? 4, 4.' )] 

xmn ilixrn i|/xn mn vm vn mn um un ymn ym jjiyn 

+2B, , (a~ l V 4' 4 +b _1 U 4 4' )(a~ 1 V 4. 4. +b' 1 Y 4, 4,' ) 

66 mn vm vn mn um un ymn 4 ym \Jj yn xmn \jixm \Ji xn 

+(D . , / a 2 ) (Y 4! 4, ) 2 +(2D 19 /ab) (y *• 4, ) (Y 4, 4' ) 

11 xmn i];xm <j/xn 12 xmn xm i|/ xn ymn ^ ym yn 


+(D 09 /b 2 ) (Y 4, 4,' ) 2 +25 1 ,(a" 1 Y 4' 4 )(a _1 Y 4' 4 

22 ymn \Jj ym \j| yn lb xmn xm | xn ymn ^ ym ^ yn 

+b" 1 Y 4. 4.' )+2D.,(b -1 Y 4 4’, )(a -1 Y $’ 4 +b‘ l 'l' 4, 4’ ) 

xmn i|/ xm \J; xn 26 ymn ij; ym tj/yn ymn ij;ym ij/ yn xmn ^ xm ^ xn 

+6,, (a 4! 4. +b ^Y 4 4! ) 2 +K A (a 4* 4 +Y 4 4 ) 

66 ymn \|fym ijj yn xmn y xm ij/xn 55 55 mn wm wn xmn ^ xm ^ xn 

+2K., -A. - (a -1 W 4’ 4 4, 4, )(b _1 W 4 4’ +Y 4 4 ) 

45 45 mn wm wn xmn djxm | xn mn wm wn ymn 4ym ilt yn 


- 1 . 


+K // A,.(b *W 4 4' +Y 4, 4 ) 2 | ab da dg 

44 44 mn wn wn ymn tyym ljlyn J 

M N i i 

-(iii/2) y y r r fm ecu 4 4 ) 2 +cv 4 4 > 2 +cw 4 4 > 2 ] 

L> L> J J L o mn um un mn vm vn mn wm wn 
m=l n=l 0 0 

+2m,[(U 4 4 ) (Y 4. 4. )+(V 4 4 ) ( w 4, 4, )] 

1 mn um un xmn i|xm tyxn mn vm vn ymn \|f ym yn 
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+ra, 


where 


SL/a(J kt 


[ (Y ft. ft, )^+(Y ft, ft. ) ]1 ab dcy dp 
xmn \|ixm \|ixn ymn tyym i|iyn J M 

,1 „1 


M N 


- 5 Y f f q W ft ft ab da 
Lj L> J J mn wm wn 


dp 


(D-l) 


m=l n=l ° ° 


ft* = fba = aft /ap 

mm n n 


D2 . Equations for Minimizing the Energy Difference 


M N 

<b/a)i u 1 1 

m=l n=l 


u I 4 


$ ft. dc* f ft ft ,dp 
„ um uk J un u-r, 

O O 


M 

N 


- \ 

v 


+A 12 L 

L 

V 

mn 

m=l 

n=l 



M 

N 

+(b/a)A 16 

I 

V 

L 


m=l 

n=l 

M 

N 



+A 16 L 

/. 

U 

mn 

m=l 

n=l 


M 

- V 

V 


+A 16 L 

L 

U 

mn 

m=l 

n=l 



M 

N 

+(a/b)A 26 


r- 

I 

> ’ 

< — i 


m=l 

n=l 

M 

N 


T 

v 


+A 66 L 

/ 

V 

mn 

m=l 

n=l 



M 

N 

+(a/b)A 66 

T 

V 
■ L 


m=l 

n=l 


M 

N 

+(b/a) ^ 

y 

V 

/ 


m=l 

n=l 


vm uk 


J 


ft t'.dff ft' ft ,dp 
— - 1 vn uf. 




® J ‘ 


$ ft . dry ft' 6 dp 
1 'um'uk J un uf. 
o o 


,1 




ft' ft .da ft ft', dp 

um uk 


un uf. 


r 1 

* • 


xmn 


ft' ft' dor ft ft ,dg 
J ijprm uk J \jfxn u-C 
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M N 


+b 10 y y v r * . da r 

12 L Li ymn tyym uk J ijiyn uf 


m=l n=l 


- ?■ v f 1 r 1 

+ < b/a)B i6 n 4 1 ZVJ o V;k d “J 

n=l c ° 


M N 


+B 


16 


y y y r* ^.d ff f *• $ .dp 

2_j 2_i xmn J \[ixm uk J ijlxn u-u 

m=l n=l 

M N i i 

+b 1a y y t r*: « , r *, #\dp 

16 Z_i 4., xmn J tarn uk J (bxn u /. 

° o ¥ 

m=l n=l 


M N j ! 

+(a/b)B 9A y J f f *, * .dof *' *'dp 

La Li ymn J tyym uk J ^yn u-f, 

m=l n=l 


+B 


66 


26 
M N 

n 

m=l n=l 


r 1 p 1 

Y $! ^ .da | i $' .dp 

ymn J q \|/ym uk i)/yn ut 


M 


+(a/b) B y yy I * *,da I dp 

66 L> ' • xmn J vjjxm uk J tyxn u-f. 


J VAA J 


m= 1 n=l 
M N 


,1 


2 v v r r 

-ab 0 / ))mU $ $ . dor $ $ .dp 

LI. o mn v um uk J un u ‘ 


m=l n=l 


M N 1 1 

2 v v r r 

-ab iXi ) ) m v 4 $ .da j 9 $ .dp = 0 

Li L 1 xmn •• tyxm uk ijjxn u/. 


m=l n=l 


(D-2) 


There is another equation, for 3L/&V ^=0, which is analogous to eq. (D-2) 
with the following substitutions: 

u*» v , U4—V, x-*y, 1*«»2, a**-b. (D-3) 


81 



M N 


~ w 1 1 

= (b/a)K A ss Y Yw f 4' *'da f 4 4 ,dB 

K’- b5 55 L-' Li mn J wra wk J wn w-f- w 


m=l n=l 


M 

+ bK 55 A 55 l 


f* 


m=l n=l 
M N 


* | 4. 4’ da f t. $ $ ,dp 

xmn J wxm wk J |xn wxn w-f. 
o o 


M N 1 i 

+ka y yw r $ 1 $ do- r § $ 

45 45 L-i 4j mn J wm wk J wn ' 


m=l n=l 


«*A* l kj 

m=l n=l 


* d0 
wn w-C, 


1 r 1 

$ 4', da 4' $ .dg 

wm wk J wn w-C. 


M N i l 

+aK,,A y y? f 4. 4 , da f 4, 4 

^5 45 L- L- xmn J ijtxm wk J \|xn ■ 


m= 1 n=l 
M N 


.. y t d$ 

\jixn w-f. 


1 1 

+b **5*45 l bynJ^ynK^l Vn». 


m=l n=l 
M N 


+<a/b) \ 4 v i k„j l 


ljiyn w-C. 


J 


d(3 


4 4 , da f 4' 4' .dR 

wm wk J wn w-f- 


m=l n=l 
M N 


o 


+aK /A/ ) y* r $ . $ , da f $ 4. 

44 44 4. ymn J ^ vjfym wk j i|iyn wf, 


m=l n=l 
M N 


'..dp 


n w i ^ 

-ab u< / ) m W 4 4 . do- | 4 4 „d6 

4. o mn J wm wk J wn w-f. 


-ab q 


..:=1 n=l 
M N , 

i n 

m=l n=l 

M N 


J V B ' 0 

o o 


(D-4) 


1 1 

sL/8 \kt ■ (b/a > s n I 


m=l n=l 
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formations 
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M N 
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(D-5) 


is an analogous expression for dL/<3Y ^ =0 , using the trans- 
(D-3) , plus the following additional one: 


45-*-*- 55 


(D-6) 
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APPENDIX E 


DERIVATION OF THICKNESS -SHEAR FACTOR FOR THE THREE-LAYER, 
SYMMETRICALLY LAMINATED CASE USING THE DYNAMIC APPROACH 


El. Dynamic Elasticity Analysis of an Individual Layer 
Undergoing Pure Thickness-Shear Motion 


For an individual layer undergoing pure thickness-shear motion 

in the xz plane, the only non-zero strain component is e xz - Now it 

is assumed that the layer is orthotropic, so that the only non-zero 

stress component is o , given by equation (49) . Then the general 

X z 

stress equations of motion for three-dimensional, dynamic elastic 
theory reduce to the following two equations: 


o 

xz ,z 


pu 


’tt 


0 

xz ,x 


pw 


tt 


(E-l) 


The longitudinal thickness-shear strain can be calculated from 
the displacements (u,w) in the x,z directions, respectively, as follows: 


e 


xz 


+ w , 


(E-2) 


The following displacement equations of motion are obtained by 
substituting equations (49) and (E-2) into equations (E-l) : 


2 . 

c s (u * 


zz + w *xz^ U * tt 


C^(u, + W , ) = w , 

s xz ’xx ’tt 


(E-3) 
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where c Is the shear wave-propagation velocity defined by: 
s 

S ' C 55 /( > <E - 4) 

However, for pure longitudinal thickness-shear motion, the dis- 
placements are independent of axial position x. Thus, all derivatives 
of displacements with respect to x vanish and equations (E-3) reduce 
to the following single expression: 



Equation (E-5) is the familiar, one-dimensional wave equation, 
which is solved easily by the separation-of -variables method, with the 
following solution for simple harmonic motion: 

u(z,t) = (A cos Q z+B sin Q z)(C, cos ht + CL sin u>t) (E-6) 

S “ i c 

where Cl s uu/c , A and B are constants depending upon the boundary 
s s 

conditions, and and are constants which depend upon the initial 
conditions. 

E2. Dynamic Elasticity Analysis for Three-Layer, 
Symmetrically Laminated Case 

Here we consider the special case of a three-layer, symmetrically 
laminated member, in which the two identical outer layers are designated 
by superscript 1 and the middle layer is denoted by superscript 2, as 
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shown in figure E-l. The proper boundary conditions for the anti- 
symmetric modes in pure thickness-shear motion are 


U (1) (z 2 ,t) = u^ (z 2 ,t) ; u.^^Zj.t) = 0 ; 
u^ 1) (z 2 ,t) = u, 2 (2) (z 2> t) ; U (2) (0,t) =0 


(E-7) 


where z^ and z 2 are dimensions shown in figure E-l. 

Substituting equation (E-6) into equations (E-7) yields the follow- 
ing set of expressions: 


A cos f/^z„ + B, sin CL^ z- = A cos z 

1 s 2 1 s 2 1 s 2 

+ B 2 sin CJ s (2) z 2 ; - sin 0 g (1) z 2 + fl g (1) Bj cos n g (1) z 2 - 


,( 2 ) 

A sin n (2) 

z„ + n (2) B 0 COS n (2) z. ; -Q (1) A, sin 

n (1) 2l 

s 

2 s 

2 s 2 s 2 s 1 

s 1 

,<» 

s 

B, cos 
1 S 

Z 1 = 0 ’ A 2 = 0 

(E-8) 


or 


cos n ( ^ z 0 
s 2 


sin z 0 
s 2 


• ->(2) 

- sm Q v 'z 
s 2 


-n (1) sin n (1) z 2 n (1) cos n (1) z -n g 2) cos n (2) z 
s s is s z s s Z 


. n <‘> or>», o‘»~. 



( \ 
A 1 


0 ' 


B 1 

II 



0 


.V 


O 

/ 


}(E-9) 


The homogeneous system of linear algebraic equations (E-9) has a 
nontrivial solution if, and only if, the determinant of its coefficient 
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matrix is equal to zero. The resulting determinantal equation has as 
its solutions the roots of the following transcendental equation: 


tan n< 2) a 2 - - n< 2) cot (n' 1 ' t 2 - 


(E-10) 


(k) (k) 

From the definitions of Q and c g , we have 


= c s (2) / c s 1} = (C 55 ) / C 55 ) )^(P (1) /P (2) ^ = (E-U) 


where 


B = c (l) /c (2) 

P " C 55 /C 55 


• » _ JD/J 2 ) 

; R = p /p 


(E-12) 


Then equation (E-10) can be expressed as follows: 


(R/p) % tan C(C 2 n < 8 1) * 1 0/R) % 3- cot [(l-C^O^ 


(E-13) 


where 


C 2 * * 2 /*! 


(E-14) 


E3. Dynamic Analysis of a Symmetrically Laminated Timoshenko 
Beam Undergoing Pure Thickness-Shear Motion 


Here we consider a symmetrically laminated Timoshenko beam . The 


A beam exhibiting both thickness-shear flexibility and rotatory 
inertia is generally referred to as a Timoshenko beam (refs. 41,42). 
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axial and thickness (or depth) directions are designated as the x and z 
axes, respectively. Such a beam could be analyzed as a special case of 
the laminated plate theory presented in Section 11 by merely deleting 
all derivatives with respect to y. However, the beam case is so much 
simpler and pure thickness-shear motion is such a simple type of motion; 
therefore, it was decided to make an exact analysis for the present case. 

The following kinematic relations hold throughout the entire 
thickness of the laminate: 


h =ili t ;e = w + \|i 

XX T X,X XZ 0,X X 


(E- 15) 


The following stress-strain relations are applicable to a typical 
layer "k": 


(k) _(k) _ p(k) \ *'/ 

o.... = E,, £.... = E, , h Z ; a - C e 


• a (k) = c (k) 

XX “11 'xx “11 "xx “ ’ xz 55 °xz 


(E-16) 


The bending moment and shear force, expressed on the basis of a unit 
width as in plate theory (Section 2.3) are: 


M 


xx 


J 


_ V 

■ L 

k=l z k-l 


(k) 

0 z 

XX 


dz 




3 (k> 

xz 


dz 


(E-17) 


k-1 k-1 


Substituting equations (E-15) and (E-16) into equation (E-17) and 
introducing the shear factor K 55 as a correction factor to be determined 
later, one obtains: 


M 

xx 



X ,x 


= K a (w + 

55 55 o,x 




(E-18) 
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where 


“ Z k 




k=l k-1 


(E-19) 


The equations of motion for a symmetrically laminated Timoshenko 
beam are identical in form to those governing a homogeneous Timoshenko 
beam (refs. 41,42), namely: 


Q = m w 
x,x o o , tt 


M - Q 

XX, X X 


m„ 


f x,tt 


(E-20) 


where m^ and m 2 are defined in equations (28) . 

Inserting equations (E-18) into equation (E-20), one obtains the 

following set of two coupled equations of motion in terms of the generalized 

displacements w and dj : 
o x 




(w + 
o.xx 


*x 



m w 
o o , tt 


(E-21) 


11 


X ,xx 


- K 55 A 55 


(w + ) = 

O ,X X 


2 w x,tt 


(E-22) 


For pure thickness-shear motion, w and ^ are independent of axial 
position x, so that equations (E-21) and (E-22) uncouple and become: 


=0 

o , tt 


(E-23) 


m. 'll _ + A__ ♦ = 0 

2 x,tt 55 55 x 


(E-24) 


Since equation (E-23) does not contain K cc and since w is not 

55 o 
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in equation (E-24) , we have no further need for equation (E-23) . 

For steady-state harmonic motion, the solution of equation (E-24) 
can be expressed as follows: 

^=? x e iUJt (E-25) 

where ^ is a constant. 

Substituting equation (E-25) into equation (E-24) , we are led to 
the following relationship: 

oi 2 = K 55 A 55 /m 2 (E-26) 

This equation is applicable to any symmetrical laminate. For the special 
case of a three-layer one, using the notation depicted in figure 10 and 
the definitions of and m^ from equations (E-19) and (28), one obtains: 


a 2 = 3K 35 z“ 2 (C^ ) /p (1) )[(C 2 /p) + 1-C 2 ] 

• [(c 2 /R) + 1 - C 2 3" 1 

where 3,R, and £ 2 are as defined previously. 


(E-27) 


E4. Determinations of the Thickness-Shear Factor 


The longitudinal thickness-shear factor, K,.,., is determined implicitly 

2 

by equating a < associated with the lowest non-trivial solution of equation 
(E-13) to that given by equation (E-27) . 
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APPENDIX F 


COMPLEX STIFFNESS COEFFICIENTS FOR A LAMINATE HAVING 
ALTERNATING PLIES OF TWO DIFFERENT COMPOSITE MATERIAL^ 

FI. Introduction 

The theory presented in Section II can be used to calculate the 
stiffness and damping coefficients for an arbitrary laminate, i.e. one 
consisting of any number of plies of any thickness, material, and orienta- 
tion. However, in many cases of practical importance, laminates are de- 
signed to have many plies of two alternating materials. In such cases, 
the approximate approach used in refs. 29-31 is sufficiently accurate. 
Although this approach was originated for application to wave propagation in 
an infinite medium, it is applicable also to plates. In the latter con- 
figuration, it has been verified experimentally in several instances. By 
comparison with static experimental results, Rose and Tshirschnitz (refer- 
ence 92) found that it gave good predictions of the in-plane elastic 
modulus and in-plane shear modulus. Also Achenbach and Zerbe (ref. 32) 
found that it gave a frequency versus wavelength relationship for longi- 
tudinal vibration of a laminated beam which was in excellent agreement 
with experimental results. 

In all of the work mentioned above, all of the individual layers were 
isotropic. However, in many cases of increasing importance, at least one 
of the sets of layers may be made of composite materials. Some examples 

*After completion of this derivation, the work of Chou et al.. (Reference 1^9) 
came to the attention of the authors. In their work, Chou et al. derived 
purely elastic equations which are analogous to the complex ones derived here . 
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are as follows: 


1. Alternating layers of an isotropic material and an orthotropic 

material. Examples: armor plate consisting of alternating layers of a 

hard ceramic material (isotropic) and a lossy glass fiber-epoxy matrix 
composite (orthotropic) ; a laminate consisting of alternate layers of high- 
modulus orthotropic composite material (such as boron-epoxy or graphite- 
epoxy) and low-modulus, high-damping polymer (to increase the damping cap- 
acity of the laminate) . 

2. Alternating layers of two different orthotropic composite materials. 
Example: Boron-epoxy (for high stiffness) and glass-epoxy (for low cost). 

The lamination scheme may be either parallel ply (unidirectional) or cross 
ply. 

3. Alternating layers of the same composite material and thickness, 
but oriented alternately at +0 and -0, where 0 < 0 < 90°. This is the so- 
called angle-ply lamination arrangement, which is very popular in a variety 
of aerospace structures. 

An original derivation is presented here which is applicable to 
determination of both the stiffness and damping of the above three classes 
of laminates. It may be considered to be a generalization of the work of 
refs. 29-31; necessarily the resulting equations reduce to theirs in the 
case when both materials are isotropic. 

F2 . Analysis 

The bases for the present analysis, as well as that of refs. 29- 
31, are the following hypotheses: 

1. Strains in the plane of the laminations are equal. 
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2. Stresses in the direction normal to lamination planes are equal. 

Hypothesis 1 leads to the Voigt upper-bound estimate (reference 93) 
of the equivalent macroscopic properties of a two-phase material having 
arbitrary geometrical configuration of individual constituents. Except 
for the presence of an additional Poisson’s ratio effects term, this is 
known also as the "rule of mixtures". Hypothesis 2 leads to the correspond- 
Reuss lower-bound estimate (reference 94), which is almost identical with 
the so-called "inverse rule of mixtures". 

Here we consider a medium consisting of repeating alternating 
layers denoted by "a" and "b", as shown in figure F-l. The plane of 
the laminations is designated as the xy plane and z is the normal to this 
plane. Then the two hypotheses mentioned above can be stated mathematically 
as follows in contracted notation: 


(a) (b) 

e. - Cj 



(j-1,2,6) (F-l) 


0 . 

1 


(a) 



( i = 3 ,4 , 5) (F-2) 


where the gj are strain components, the are stress components, and 
superscripts (a) and (b) refer to layers "a" and "b". Subscripts 1 and 2 
refer to normal strain (or normal stress) in the plane of the layers, 3 
refers to thickness-normal effects, 4 and 5 refer to thickness shear, and 
6 refers to in-plane shear. This notation is consistent with that most 
widely used in the field of composite-material mechanics (ref. 4) but 
differs from the mixed notation used in the body of this report. 

As a consequence of the complex, orthotropic version of Hooke's law 
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and static equilibrium, equations (F-l) and (F-2) imply the following re- 
lations : 


o. 

J 


H W o. (a) + H (b) 0. (W 
J J 


( j = l »2 ,6) 


e 


i 




„(b) (b) 

h e . 


(i-3,4,5) 


(F-3) 


(F-4) 


where and are the thickness fractions of layers "a" and "b", 

respectively. Thus, 


,<a) 


+ H< b) - 1 


(F-5) 


The complex, orthotropic version of Hooke's law, alluded to above, 
holds for each type of layer as follows: 


>> 


" p(k) 

Si 

M 

°12 

s(k) 

13 

0 

0 

p(k)' 

C 16 


f \ 
e i 

(k) 

°2 


p(k) 

°12 

p(k) 

22 

p(k) 

23 

0 

0 

p(k) 

26 


e 2 

°3 

< 

>• 

s( k) 
C 13 

? (k) 

L 23 

p(k) 

C 33 

0 

0 

p(k) 

36 


(k) 

e 3 

> 

°4 


0 

0 

0 

p(k) 

C 44 

p(k) 

C 45 

0 


e (k) 

e 4 

°5 


0 

0 

0 

p(k) 

°45 

p(k) 

C 55 

0 


(k) 

e 5 

0< k > 
l 6 > 


p(k) 

C 16 

p(k) 

26 

? (k) 

C 36 

0 

0 

p(k) 

C 66 


e, 

6 J 


(k=a,b) 


(F-6) 


-(k) 

where the C. . 

J-J 


are the complex stiffness coefficients. 
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S' 


v 3 J 


a il 

6 i2 

r 1 

L 13 

c 1 1 
13 

c 1 
C 16 

a { 2 

^22 

c 1 

23 

C 1 1 
23 

c 1 

°26 

p(a) 

C 13 

p(a) 

C 23 

p(a) 

°33 

0 

p(a) 

36 

p(b) 

C 13 

~<b) 

23 

0 

s(b) 

C 33 

? (b) 

C 36 


"1 

r«x ' 


e 2 

< 

(a) 

e 3 


(b) 

e 3 

- 

6, 

^ 6 J 


(F-7) 


where 


c h 


C! . 
ij 


C! ! 
ij 


H <a> C U) + H (b) £? b) ; i,j 
1J 1J 

H (a) C^ a) ; i,j = 13, 23 

H (b) 5 ^(b). ... . 13> 23 


= 11, 12, 16, 22, 26 


Combining equations (F-4) and the last two of equations (F-7), we 
get the following relations: 


(a) r_(b) H (b), = (a) nV u / /pV°/ p \ \ 

6^ ”” 1^23 ^3 — ^ '^13 * ^*13 '^3 ^ i ^ 


(b) ,^(a) p (b) . 

23 " 23 ' 6 2 


H <W (C ( 3 b) 6 - C< a) ), 6 ] [H <a) C<»+ H (b) c'f]' 1 


(F-8) 


(b) _ rM 


.(a) .(b) -(a) 


1^33 ^ a / 6 1 " v'-'o'l ™ U,* ' 


(a), -(b) -(a). 


'13 13 /e l 


23 23 fc 2 


■H (a> (C< b) - c‘ a) ) C6 ] [H (a >c' b> + H (b > c'Jr 1 


(F-9) 


Substituting equations (F-8) and (F-9) into the first of equations 
(F-7), one is led to the following result: 


°1 = C il e l + C 12 e 2 + C 13 e 3 + C 16 e 6 


(F-10) 
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where 


p _ H (a) p(a) (b)-(b) £ ,p(a) = (b ).2 

C 11 ^ H C 11 + H C U “ K C (C 13 ' C 13 } 


(F-ll) 


C i2 - H (a) C^ ) + » (b) C< b) - K c (C< a) - C^HC^-cj^) (F-12) 


5 13 " t5i3 ) < H<a)/e 33 )) + C^(H (b) /C<W)l L. 


(F-13) 


C = p( a ) t u(k)r(b) V “(b)v - .v 

C 16 H C 16 +H C 16 " K c (C 13 C 13 )(C 36 °36 > (F ~ 14) 


K c - [(C^ ) /H (a) )+(C^3 ) /H (b) )]- 1 ]" 1 


(F-15) 


L c = L(H (a) /C^ ) ) + (H (b) /c^ ) ) J ' 1 


(F-16) 


Analogous equations to equations (F-10) - (F-14) can be obtained 
readily by interchanging the roles of subscripts 1 and 2 . 

From the last two of equations (F-7) , one gets: 


9rr _ /-(a), p( b )\ _ . r p(a) + + p( a ) (a) + ?( b ) (b) 

2ct 3 (C 13 + C 13 )e l + °23 C 23 * 2 + C 33 e 3 °33 e 3 


+ (£ 36 } + a 36 )} e 6 


(F-17) 


Then substituting equations (F- 8 ) and (F-9) into equation (F-17), one 
finds the following result: 

°3 = ^13 e l + ^23 e 2 + ^33 e 3 + ^36 e 6 (F-18) 

where C 33 and ^3 are given by equation (F-13) and an analogous equation, 
and 



= L 

c 


(F- 19) 
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(F-20) 


£ 36 ' C'W 


(H (a) /C< a) 


) + C<W(H <b >/C<W) ] 


Still another set of equations analogous to equations (F-10)-(F-14) 
can be obtained by interchanging the roles of subscripts 1 and 6. 

As the first step toward determining expressions involving the com- 
posite coefficients with subscripts 44, 45, and 55, it is expedient to in 
vert the pertinent equations in set (F-6) , thus arriving at the following 
expressions: 



(F-21) 


where 


5<k) = r (k) /A (k) 

S 44 ~ C 55 /A 


S (k) = C /A (k) 
s 55 C 44 /A 


-(k) = x(k),-(k) 

45 - " 45 /A ’ 


rOO . = 00 c (k) ffi< k V 

A - C 44 C 55 (C 45 ) 


(F-22) 


Inserting equations (F-21) into equations (F-4) , we arrive at the 
following expressions: 


where 



(F-23) 



H (a) S<f + 
44 


H (b) S^ 

44 


(F-24) 
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There are three expressions, analogous to equation (F-24) , which are 
obtained by merely substituting subscripts 45 and 55, respectively, 
for subscript 44 in equation (F-24) . 

We now have derived expressions for all twenty of the laminate 
coefficients which appear in equation (F-23) , above, and in the follow- 
ing expression: 



C 

C 

C 

C 


11 

5 12 

5 13 

S6 

12 

^ 22 

^23 

a 26 

13 

5 23 

£ 33 

a 36 

16 

£ 26 

5 36 

^66 


(F-25) 


The effective density of the laminate is given by the following 
expression: 


p 


H< a) p (a) + H< b V b > 


(F-26) 
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APPENDIX G 


IDENTIFICATION OF INTEGRAL FORMS 
Gl. Trigonometric Integrals 

The integral forms for three boundary conditions are tabulated in 
the following form, where m,k,n,-t. = 1 to 2: 



Integral 

Evaluation for 

Boundary Condition 

Listed 


Form 

Simply Supported 

Clamped Edges 

Free Edges 

I, , 

= $ $ dor 

% ; 

m=k/0 

k ; m=k/0 

o 

ii 

B 

CM 

\ 

ro 

lmk 

j um uk 
o 

o ; 

m^k 

0 ; m^k 

1 ; m^k 

hnt 

=? f i $ ,dp 

J un ut 
o 

% ; 
o ; 

n=Lfk) 
n tL 

h ; n~li*0 
0 ; n^L 

3/2 5 n«t*> 
1 ; ntl 

I- , 

= $ § . da 

\ ; 

m=k^0 

* lmk 

X lmk 

2mk 

J vm vk 
o 

o ; 

m^k 

1 2nt 

- f $ $ ,d0 

J vn v-f- 
o 

% ; 
0 ; 

n =t#) 
n 

1 ln-t 

X lnt 

*3mk 

= $ $ . da 

J um vk 
o 

0 


I lmk 

* lmk 

I 3n-t 

= $ $ dp 

.1 un v-t 
o 

0 


^Int 

1 ln-t 


= [ $ *’ dor 

•J um uk 

kn/2 

; m=k/0 

0 ; m^k^O 

0 ; m=k#0 

J 4mk 

0 ; 



o 

m^k 

0 ; m/k 

0 ; m^k 


r 1 

nn/2 

; n=£?*0 

0 j n sa L?^0 

0 ; n=L ^ 0 

*4nL 

S $ #' dp 

J un v-C 
o 

0 ; 

n^-t 

0 ; nM 

0 ; nH 
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I 


I 


X 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


Integra 1 

Evaluation for 

Boundary Condition 

Listed 


Form 

c . s [* $ $' .do 

5mk J vm uk 

o 

Simply Supported 

Clamped Edges 

Free Edges 


krr/2 ; m=k^0 
0 ; m^k 

*4mk 

*4mk 


«. , » f « *’dp 

5nt 0 vn ut 

o 

nrr/2 ; n=t^0 
0 ; n^t 

X 4nt 

I 4nt 


= f 6 $ 1 do 
6mk J um uk 
o 

0 

^4mk 

^ 4mk 


6n X s J *„n‘^ d e 

o 

0 

*4nt 

X 4nt 


f l 





= f #’ dor 

7mk J vm vk 

o 

0 

^4mk 

^4mk 


r l 





7 . s * *' dp 

7n{. J vn v-L 

o 

0 

T 4nt 

I 4nt 


r 1 

mkn^/2; m=k^0 

2mktT^ ;m=k^0 

2mkrr^ ; m 

i=k^0 

a . S $’ a’ do 

8mk J um uk 

o 

0 ; m^k 

0 ; m^k 

0 ; m^k 


8nt S I *un $ ut dP 
o 

2 

ntrr / 2 ; n=t^0 
0 ; n^t 

2 

2ntrr ;n=t^0 
0 ; n^t 

2n-t-TT^ ; 
0 ; n^t 


= r $* $• do 

9mk J vm vk 
o 

mkrt^/2 ; ra=k^0 
0 ; m^k 

^8mk 

^8mk 


. 5 f f $’.dp 
9nt J vn vt 
o 

Mn2/2 ; n=t?*0 
0 ; n^t 

I 8nt 

I 8n^. 


r f $ 1 $ 1 dO 

10mk J um vk 
o 

0 

^8mk 

^8mk 


,1 





in , s *’ 

lOnt J un vt 

o 

0 

r 8nt 

X 8nt 


f 1 

~ I $ 1 $ do 

llmk J o t | ixm wk 

0 

See G2 for 
c lamped case 

0 


. = f $ $ , dB 

lint J tocn wL 
o 

Jg ; n=t/0 
0 ; n^t 

See G2 

0 
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Integral 

For© 


$ dor 


, rnn dj ri on List ed 

Simply Supported CUrnP 

% ; m-kW I llmk 

llmk 

n : ©^ k 


ktr /2 ; 


See G2 


See G2 


_ r$< $ /P 0 
' J 0 ^* n Wt 

- [ l $ *\ 6 * ° 
h5mk"J 0 ^ Wk 

nfl/2 *, n-t^° 

I,e n ; “ ! ^liiyn^w-t ^ 0 ; 


1 13mk 


i »i a, dor 
wk 


° nn/2 ; n=/ ^° 

- I 0 ; n*. 


hen-t ^ J 0 w n 

_ f 1 * $ . dof 

Il7n.1t " J- W ® Wk 


• m=k ^0 


k ; ne -^° 


I 17nt * I 0 ^^ d6 0 ; ^ 

_ fv $ do? 0 

1 lRn.lt " J ^ 


Hsnt ' 


$ ., d P 


i 



Integral 

Form 


Evaluation for Boundary Condition Listed 
Simply Supported Clamped Edges Free Edges 


19mk 


s f *' i.do 0 

J wm wk 
o 

I 1Q , s f $' $ f dg 0 
19n-u J wn w( ^ 


18mk 


18n.t 


20nf 


21mk 


21nt 


22mk 


22n-f, 


23mk 


23n-f. 


24mk 


25mk 


2 3n-t, 


18mk 


18n-t, 


*1 

mkn^ / 2 ; m=k^0 


2 

mkn /2; m=k?k) 

I f'da 


C Z (C Z -2) 

J wm wk 
o 

0 ; m#k 

mm mm 

0 ; m^k 

r 1 

ntn^/l ;n=t,^0 


n£n2/2 ; n^^O 

$' $',dp 
J wn wl 
o 

0 ; n^t 

C Z (C Z -2) 
n n n n 

0 ; n^t 

r 1 

$ , $ , , da 

\ ; m=k^0 

I, , 

T 

fxm \)jxk 

0 ; m^k 

lmk 

lmk 


k ; n=f,^0 
0 ; n^i 

T lnA 

I lni 

f 1 

\ 5 m=k#) 



J o ^>[(ym^\Jiyk da 

0 ; m?*k 

^lmk 

^ lmk 

r 1 

^ ; n=-t,j*0 



J^tyyn^y^ ^ 

0 ; nU 

^n-f. 

r int 

r 1 




$ , $ . , da 

0 

I , 

I 

J fxm \[iyk 

lmk 

lmk 

r 1 




J ^ilfxn%y-C, d ^ 

o 

0 

I ln-t 

I ln-f- 

p* 

kn/2 ; m=k^0 



J o ^ij/xm%yk dc * 

0 ; m^k 

*4mk 

*4mk 

r 1 

nn/2 ; n=£^0 



$. .dp 

J \jtxn ify’C 

0 ; nM 

*4nt 

^4n-t 

pi 

kn/2 ; m=k^0 




0 ; m^k 

*4rak 

*4mk 

r 1 

nn/2 ; n^^O 



J o ^yn%x^ dP 

0 ; n 

*4n-f. 

^4n-t 
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Integral Evaluation for Boundary Condition Listed 


Form 

Simply Supported 

Clamped Edges 

Free Edges 

X 26mk ~ 

f 4. 4! , doc 

.1 ijtxm \Jjxk 

0 

*4mk 

*4mk 


r 1 




I 26nt 

l o V»V e 

0 

■*■ 40 /. 

I 4n-f. 

I 27mk 

I o $ 'liym^i(fyk d0 ' 

0 

*4mk 

J 4mk 


A 




I 27n-t ~ 

J 0 Sy»%y* dp 

0 

^4n-4 

X 4nt 


r 1 

2 

rnkn /2 ; m=k^0 



1 28mk s 

$ ’ 4 ' da 

ijixm iyxk 

0 ; m^k 

*8mk 

^8mk 


r 1 

ntrr^/2 ; 



I 28nt ~ 

J o ^xn%xJ d ^ 

0 ; 

I 8n-f, 

I 8ni 


r 1 

mkn^/2; m=k? t 0 



^29mk 

I 4! 4 ' da 

J \|iym ipyk 

0 ; m^k 

^8mk 

*8mk 


,1 

ntn^ / 2 ; n=£^0 



X 29n-t = 

4’ 4! ,dp 

•' 0 ^yn tyy*- 

0 ; nH 

I 8n-(. 

l 8nl 


f i 




I 30mk = 

4 ! 4' .da 

J \jixtn \|iy k 

0 

*8mk 

I 8mk 


r 1 




X 30nt = 

$' 4' dB 

J 9 ^xn ijiyf P 

0 

X 8n/. 

l 8nl 


r 1 

\ ; m=k^0 

^ lmk 

1 lmk 

^lmk 

1 o ^\)ixm^uk da 

0 ; m^k 




pi 

\ ; n-*,#> 



X 31nf “ 

1 4 4 p dB 

•j yxn ut 

0 ; nil 

^ln<t 

J ln4. 


f 1 

\ ; m=k#D 



X 32mk ~ 

4 4 , . da 

.1 um ijixk 

0 ; mj^k 

*lmk 

I lmk 


r 1 

k ; n=t^0 



X 32nt 5 

4 4 .dp 

J ^ un ij/xt. 

0 ; n£f. 

1 ln-C 

1 ln-t 
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Integral Evaluation for Boundary Condition Listed 


Form 

Simply Supported 

Clamped Edges 

Free Edges 

f 

,1 




X 33mk i 

$ , dcr 
tyxm uk 

0 

X 4mk 

X 4mk 

1 33n-t " J 

.1 

§' $ „df} 

Axn ul 
o 

0 

X 4n-{. 

X 4n-f, 

f 

.1 




X 34mk J 

$ 1' do 

pm uk 
o 

0 

X 4mk 

X 4mk 

f 

.1 




I 34nt = J 

i. *'.d0 
yxn u c 

0 

X 4n£ 

X 4n-t 

| 

d 




I 35mk = J 

$ . ^ .do 

. tyym uk 

0 

X lmk 

X lmk 

I 35nl s J 

r 1 

W 6 

o 

0 

I ln-6 

1 ln-f. 

| 

.1 

kn/2 ; m=k^0 



X 36mk 

$ . do 
tyym uk 

0 ; m/k 

X 4mk 

X 4mk 


a 

nn/2 ; n=t^0 



I 36nl " J 

$’ $ # dp 

vlryn u-L 
o 

c 

o 

X 4n{, 

X 4n-t 


r 1 

kn/2 ; m=k/0 



X 37mk 5 * 

$ $’ da 

1 mym uk 
o 

0 ; m^k 

X 4mk 

X 4mk 


1 

nn/2 ; n*A^0 



T 37n-f. ® . 

f o W e 

0 ; n^-l 

X 4n/, 

X 4nt 


A 




X 38mk 

$ . $ . do 

! uixm vk 
o 

0 

X lmk 

X lmk 

X 38nt 5 • 

f $. $ ,dg 

J u«n vl 
o 

0 

1 ln-f. 

I lnt 


r l 

kn/2 ; m=k#) 



X 39mk 

1 $ ' $ . dor 

J ikxm vk 
o 

0 ; m^k 

X 4mk 

X 4mk 
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Integral 

Form 


Evaluation for Boundary Condition Listed 
Simply Supported Clamped Edges Free Edges 


39n l 

c 

f *; $ .dp 

1 tfixn vf 
o T 

nrr/2 ; n=£^0 
0 ; n^-t 

I 4nt 

^4n-t 

"40mk 

t 

f $ . $ do- 

1 wxm vk 
o 

krr/2 ; m=k#0 
0 ; m^k 

^4mk 

^4mk 

40n-t 

t 

A 

*. t'.dp 
fan vl- 
o T 

nrr/2 ; n=l^0 
0 ; 

X 4nt 

: 4n-t 

41mk 

= 

$ . $ , da 

1 rjrytn vk 

\ ; m=k^0 
0 ; m^k 

^lmk 

^lmk 

41 nl 


p l 

$ § ,d(3 

\|iyn v-f/ 

\ ; n=<^0 
0 ; n=/t 

1 lnt 

1 ln-C, 



,1 




42mk 

= | 

$ ' f da 
ij/ym vk 

0 

^4mk 

^4mk 

42nl. 

t 

,1 

$' $ .dp 

o tyyn v-?. 

0 

*4n-f. 


43mk 

t 

r 1 

$ $' da 

r|rym vk 

0 

^4mk 

' t 4mk 



,1 




43n-t 

i 

1 $. *' dp 

r|ryn v-L 

0 

^4n-t 

I 4n-t, 

44mk 

■J 

r 1 

$ ' $ ' , da 

fam uk 
o 

mkn^/2; m=k/0 
0 ; m^k 

*8mk 

^8mk 

44n-f. 

■J 

.1 

$! 

fan ul 
o 

n-f-Tr2 /2 ; n=f ^0 

0 ; n/l. 

I 8n-f. 

^8n£ 

55mk 

-J 

.1 

$' da 

\Jiym uk 

0 

*8mk 

'''Smk 



J 




55n-f. 


*; t'.dp 

tjryn u-c 

0 

L 8n/. 

I 8nt 



J 




56mk 

* J 

$ ' $' da 

l tyxm vk 

0 

^8mk 

^8mk 
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Integral Evaluation for Boundary Condition Listed 


Form 

Simply Supported 

Clamped Edges 

Free Edges 

56nt 1 

•1 

$! $',dp 

\lixn vt 
o T 

0 

I 8nt 

I 8nt 

* f 

57mk J 

.1 

do 

(Jfym vk 

mkTT 2 /2 ; m=k^0 
0 ; m^k 

*8mk 

^8mk 

57nt = J 

.1 

f ! §',dg 

Q fyn vt 

ntn^ / 2 ; n=t^0 
0 ; nH 

^8nt 

X 8nt 


G2 . Combination Trigonometric-Beam Type Integrals 

These integrals are related to I llmk >I Unr I 13mk > and 1^ for 
clamped and free edges of the plate mentioned in Appendix Gl. Therefore, 
these integrals were evaluated and are listed below: 

All Edges Clamped . 

r 1 

^llmk J ^tyxm^wk^ 1 * 

o 


(Z k ) 2 + (2om) 2 


|z k sinh Z^ sin2mn-2mn cosh Z^ cos 2mn + 2mn| 


fcos[2mn + Zk]-1 cos[2mTT - Z k ]-1 
t 2[2mTT + Z. ] + 2[2mTT - Z ] 


(Z ) 2 +(2mTT) 2 1 k 

k 


|z, cosh Z, sin 2mn - 2nm sinh Z, cos 2mnJ 


+ C 


, sin[2mTT-Z, 1 

{ *LL 

k l 2r2mT7 - Z ] 
k 


sinT 2mrr + Z ] , 

' L_ l 

2f2mTT + Z^] J 


I is analogous to I,. , with the following substitutions: 

lint 1 lmk 
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m ■*->- n , 


k t , a *-*■ 3 


I 


1 3mk 


r 1 2 

*. J’.do - K 

wxm wk 

o (Z, ) + (2mrr) 


0 « i Z, cosh Z, sin 2mrr 

2 . , n v 2 l k k 


- 2mrr sinh Z, cos 2raTr r + Z 
k 


, 8 Inf 2 mn - 21,1 ■!«[>"+ V 

/ k \ 2[2mTT - Z k ] 2[2mTT + Z ] 


Z k C k 


— 9 1 Z i, sinh Z sin 2mn - 2mn cosh Z cos 2mn + 2mTT| 
(Z k ) Z + (2mTT) 




,cos[2mn + Z k ]-1 cos[2ran - z k ]-l 

Z k ^k i 2 f 2nm + Z ] 2[2mTT - Z ] J 


Il3n> * S ana l°8 ous to Ij^ mk with the following substitutions: 


m n , k f , cy ► 3 
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APPENDIX H 


EVALUATION OF EXPERIMENTAL METHODS USED TO DETERMINE MATERIAL 
DAMPING IN COMPOSITE MATERIALS 


HI. Specimens and Experimental Techniques 


The following experimental techniques and specimen types have been 
used to determine material damping by experimental means (references 
48, 79) : 

1. Decay of free vibration: 

a. Torsional pendulum 

k 

b. Axially vibrating bar (reference 95) 

c . Beam 

k 

(1) Free-free (reference 96) 

k 

(2) Cantilever (references 97, 98) 

2. Resonant response: 


a . 

Torsional 

(reference 99) 

b. 

Axia 1 


c . 

Beam 



(1) Free- 

k 

free (references 99, 100) 


•k 

(2) Cantilever (references 101-104) 


Apparently only these methods have been applied to measurement of damping 
in composite materials; each is discussed briefly in succeeding paragraphs. 
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d . Cubic block 


e. Plate (reference 13) 

ic 

f. Shell (reference 105) 

* 

3. Rotating-beam deflection (reference 106) 

4. Composite oscillator 

5. Ultrasonic methods 

6. Thermal methods 

Pottinger (reference 95) measured the temporal decay of axial vibrations 
of bars (Type lb, above) of glass fiber-epoxy and boron fiber-aluminum com- 
posites in the frequency range of 1 kHz to 100 kHz. 

Among the earliest uses of free vibration to determine the logarithmic 
decrement of a composite material were the free-free beam experiments re- 
ported by Bert et al. (reference 96). The specimens were sandwich beams with 
glass-epoxy facings and hexagonal-cell honeycomb cores of either aluminum or 
glass-phenolic. Good correlation was obtained between the measured values 
of 6 and those predicted by an energy analysis of a Timoshenko beam, which 
includes thickness-shear flexibility and rotatory inertia. However, the 
complexity of the analysis prevents explicit determination of the 6's of the 
facings and core. 

Later Schultz and Tsai (references 97, 98) used free flexural vibration 
of very thin cantilever beams to determine 6 for glass-epoxy in two con- 
figurations: (1) unidirectional with the fibers at an arbitrary orientation 

0 with respect to the longitudinal axis of the beam, and (2) symmetrically 
laminated quasi-isotropic layups (layer orientations of 0°, -60°, 60°, 

-60°, 0° and of 0°, 90°, 45°, -45°, -45°, 45°, 90°, 0°). In the 9 t 0°, 

90° case, the coupling between bending and twisting significantly in- 
validated the simple flexural theory used in the data reduction. 
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Adams et al. (reference 99) used torsional and flexural resonant 
response to determine the damping of glass fiber-epoxy and carbon 
fiber-epoxy. The magnetic-coil driver was located at the one end 
of the test bar to excite torsion, Type 2a, and at the center to 
excite flexure, Type 2c(l) . Similar tests were conducted by Wells 
et al. (reference 100) on graphite-epoxy. 

Cantilever beams, Type 2c(2), were used by Keer and Lazan (reference 
101), Gustafson et al. (reference 102), Mazza et al. (reference 103), 
and Tauchert and Moon (reference 104). It is noted that Keer and Lazan 
used sandwich-type beams with glass reinforced plastic facings and core, 
and determined the resonant energy dissipated per cycle as the measure 
of damping, and found it to be proportional approximately to the square 
of the stress amplitude, i.e. Kimball-Lovell behavior (see Appendix B) . 

Clary (ref. 13) conducted resonant response tests on long free-edge 
plates (Type 2e) made of unidirectional boron fiber-epoxy. He used 
the modified Kennedy-Pancu method of data reduction, and studied the 
effect of fiber orientation. The analytical results of the present 
study are compared with his experimental results in Section 3.3. 

Bert and Ray (reference 105) carried out resonant response measure- 
ments on a free-edge, circular, truncated conical shell (Type 2f) with 
an aluminum honeycomb core and glass fiber-epoxy facings. They used 
the original Kennedy-Pancu data reduction method. 

Richter (reference 106) used the rotating-beam deflection technique 
(Type 3) to determine the damping characteristics of glass fiber-epoxy 
at low frequency (0.01 to 1.67 Hz). 

It is apparent that the most popular techniques that have been used 
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in connection with measurement of composite-material damping are the 
free-vibration decay and resonant-response methods. Differences in re- 
sults obtained by these two methods have been studied analytically by Parke 
(reference 73), using a standard linear solid, and by Heller and Nederveen 
(reference 107), using a generalized Maxwell model. 

H2. Excitation and Data-Reduction Techniques for Modal Resonant 
Response of Complicated Structures 

The experimental and data-reduction techniques for determining damping 
from free-vibration decay are rather straight forward. In contrast, a 
great variety of techniques are in use for modal response at resonance, 
especially in regard to data reduction. 

One approach is to use multiple shakers located at various locations 
on the structure and having controlled force amplitudes and phase re- 
lationships. This approach was origated by Lewis and Wrisley (reference 
108) and has been used and extended by Fraeijs de Veubeke (reference 109) , 
Asher (references 1 10 , 1 1 l),Tra i 1 1-Nash et al (references 1 12 , 1 13) , Hawkins 
et al (references 114,113), and de Vries (reference 116). The major dis- 
advantage of this approach is the complexity and expense, in time and cost, 
of installing multiple shakers. 

Another approach is to use a single excitation, even when the structure 
is complicated. In conjunction with this method, there are a variety of 
resonance-determination methods; these are discussed in the paragraphs which 
follow. Comparative critical evaluations of a number of these methods 
have been presented by Bishop and Gladwell (reference 117) , Pendered and 
Bishop (reference 41) , and Pendered (reference 118) . 
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Peak-Amplitude Method . - In this simple, classical method, resonance 


is defined very simply to occur when the amplitude response reaches a 
peak. The major disadvantage of this method is that it is limited to 
lightly damped systems without closely spaced frequencies. The reason for 
this is that both of two closely-spaced-frequency modes may contribute to 
the amplitude at an intermediate frequency resulting in an indication of 
only one mode (a pseudo mode, of course). In an analytical study of a 
typical system with two natural frequencies, closely spaced, Turner 
(reference 119) showed that the peak-amplitude method results in a 1307., 
error in determination of the resonant amplitude. 

In general, there is no direct means of determining the modal damp- 
ing, i.e. the magnitude of damping associated with various vibrational 
modes. However, it can be determined from the response at the anti- 
resonant frequencies, using the equations derived by Brann et al (refer- 
ence 120). Also, for the special case of a base-excited system, the ratio 
of the tip-to-base amplitude to the base amplitude is related to the 
damping, as was shown in references 121,122. 

Peak-Quadrature-Component Method . - This method was originated by 
Stahle and Forlifer (reference 123) in 1958. In this method, resonance is 
defined to occur when the quadrature component (90° out of phase with the 
excitation) of the response peaks in either a positive or negative dir- 
ection. Damping is determined from the in-phase-component peaks which 
occur at a frequency below resonance and another frequency above re- 
sonance. Applications of the method are discussed in references 124-126. 

Quadrature-Response Method . - In this method, resonance is defined to 
occur when the in-phase component of the response vanishes. The errors 
associated with this method have been discussed by Pendered (ref. 118). 



Kennedy- Pancu Method and Its Modifications . 


This method was origin- 


ated by Kennedy and Pancu (reference 87) in 1947. Its basis is the geo- 
metry of the theoretical Argand diagram for a single-degree-of-freedom 
damped system, which is a circle (see figure H-l). The data-reduction 
procedure is to fair a circle through an Argand diagram obtained from 
experimental data corresponding to a range of frequencies. The tacit 
assumptions inherent in the method are: 

1. The system is a lightly damped, linear system. 

2. In the vicinity of a particular resonance, the extraneous con- 
tributions from the off-resonant modes are either negligible or invariant 
with respect to frequency. 

3. There is no damping between the normal modes, i.e. the inertia, 
damping and stiffness matrices are all diagonal. 

Advantages of the method are: 

1. It permits resonant frequencies, modal amplitudes, and modal 
damping values all to be determined from the Argand plot. 

2. It has been demonstrated on a comparative basis that the KP 
(Kennedy-Pancu) has better capability in separating closely spaced re- 
sonant frequencies. 

In the KP method, two ways to determine the resonant frequencies are 
possible. One way is to use the frequency associated with the point (on 
the faired-in circle) which is most distant from the real axis. The other 
way is to use the maximum frequency spacing technique, in which resonance 
is defined to occur when the rate of change of the arc length (of the 
Argand plot) with respect to frequency attains a local maximum. 

In the original Kennedy-Pancu method, the damping ratio was determined 
from the half-power points, which are two ends of the diameter (of the 
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Argand circle) which is parallel to the real axis. See equation (B-56) . 

Pendered and Bishop (reference 88) introduced an improved method, based 
upon the change in phase with respect to frequency, evaluated at re- 
sonance. See equation (B-69) . 

Gladwell (reference 127) proposed a method of determining the true 
response peaks and their associated damping ratios, using three points 
on the frequency-response curve. However, Pendered and Bishop (ref. 88) 
showed that variations of the two side-points caused changes of 135% 
in the damping coefficient of the system that they investigated. 

Keller (reference 128). Smart (reference 129) , and Pallett (reference 
130), discussed improvements in instrumentation used in conjunction with 
the KP method. In his analytical investigation of a typical two-degree- 
of freedom system with closely spaced natural frequencies. Turner (ref. 

119) indicated a 267„ error in resonant amplitude determination, which is 
a considerable improvement over the 130% error for the peak-amplitude 
method. The modal-shape-determination aspects of the KP method were 
discussed in reference 131. 

The KP method has been applied successfully to a great variety of 
vibrating systems, including: 

1. Aircraft ground and flight testing (ref. 87) 

2. Beam-type structures (ref. 88) 

3. Two-degree-of-freedom system (ref. 120) 

4. Liquid-propellant launch-vehicle axial vibration (refs. 132,133) 

5. Unstiffened stringer-stiffened, and ring-stiffened sandwich cylin- 
drical shells( refs. 103,130,134,135). 

6. Sandwich conical shell with composite-material facings (ref s . 102-104) 
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7. Composite-material plates (ref. 13). 

The last three categories listed above are particularly useful ap- 
plications of the mode-separating capability of the KP method, since they 
often have numerous closely spaced natural frequencies. Tendered and Bishop 
(reference 136) used the KP method in connection with determination of the 
dynamic characteristics of a sub-system from resonance test results on 
a complete system consisting of two sub-systems. 

Woodcock (ref erence 137) extended the KP method by providing for any 
arbitrary amount of damping of the classical viscous type and by 
making no assumptions about the form of the damping, stiffness, and in- 
ertia matrices. However, as pointed out by Nissim (reference 138), it 
is necessary that the inertia matrix be known a priori in order to deter- 
mine experimentally the complete dynamic description of the system. (It 
is noted here that the inertia matrix could be determined experimentally, 
at the cost of additional tests, by the displaced-frequency method, refer- 
ence 139) . Nissim went on to present a set of data-reduction equations 
with which the disadvantages of ref. 138 are eliminated. It should be 
pointed out that apparently no actual applications of the methods of 
refs. 137 and 138 have been documented so far. 

H3. New Equations for Applying Two Different Versions of the Kennedy- 
Pancu Method to Two More Realistic Classes of Solids 

From linear vibration theory, it is well known that the largest value 
of the damping ratio Q for which free vibrations can occur is £=£ *1; 

at higher values of £ , the motion is nonperiodic. For a standard linear 
solid, Zener (reference 140) showed that the critical loss tangent is 53°, 
which corresponds to approximately 1.183 for g . (This compares with a 
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value of 2.00 for g cr of a Kelvin-Voigt or Kimba 1 1-Lovel 1 solid, since 
equation (B-65) is assumed to hold for such materials). Parks (ref. 

73 showed the error resulting from the use of equation (B-75) for a 
standard linear solid. 

As mentioned previously in the preceding section, in applying 
the Kennedy- Pancu method, it is tacitly assumed that the inertia, damping 
and stiffness matrices are all diagonal. Although this assumption may be 
satisfied in the case of a simple structure made of one material, it is 
certainly not justified in the case of a laminated composite-material 
structure. However, if the damping is sufficiently small that it does 
not affect the modal shapes (although, of course, it affects the modal 
amplitudes) and the intermodal damping coupling is negligible, then the 
KP method should still be sufficiently accurate to justify its use 
for engineering purposes. 

It was mentioned in Section B3 that, in order for two half-power 
points to exist for a Kimball-Lovell material, the damping ratio Q 
must be less than/2/2. However, when the damping is greater than this 
value, one half-power point (denoted as 012 ) still exists. Thus, the 
method can be modified to accommodate this by using equation (B-55) re- 
written as follows: 

c = (1/2) [ (&' /it ) 2 - 1 ] ' (H-l) 

l n 

Just as the modification of the original KP method as described above, 
equation (H-l), permits its use regardless of the amount of damping, other 
appropriate equations can be derived to extend the use of the method to 
linear systems composed of materials represented by other damping models. 
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Here we present original equations for two cases: the standard linear 

solid (or first-degree hereditary mode 1) and the Biot model. Since 
these models contain two parameters rather than one as in the cases of 
viscous and Kimball-Lovell damping, more experimental data will be 
necessary, of course. 


Standard Linear Solid (or First-Degree Hereditary Model) . - 
Volterra (ref. 83) presented an analysis of the sinusoidally forced 
vibration of a single-degree-of-freedom system exhibiting first-de- 
gree hereditary damping, which is equivalent to a standard linear 
solid. The equation of motion can be written in the present notation 
as follows: 

mu + ku + h f e ^l^ [du(-r) /dr]dT = F sin out (H-2) 

J o 

o 

or 


mu 4<*^mu +(k+h)u +<x^ku = F q (uu cos out + o^sin uut) 


(H-3) 


The steady-state-response solution of equation (H-3) is of the 

form: 


u = u sin (ait - m) 
o r 


(H-4) 


where 


x /r 2 2 , 2 u 2 2 2 4 -|2 r „2 ,2 11/2 

u = (F /m) I La uu - (a. -H -jj )x -ou J +[a,H uuj r 
ootnin ij 

.[ 2 , 2 2.2 2 2 2 2 2 1 -1 

[a i (uo^ - jj ) + uu (uu -a» n -H ) J (H-5) 

f „ 2 r 2 2 , 2 2 2 . 2 I 

cp = arctan |o^H uu “ (or^-H -uu^lou -uu J j- (H 6) 


118 



2 2 
H ~ h/m , hi h k/ro 
n 


Equation (H-5) can be written in dimensionless form as follows: 


MF « [[a 2 - (o 2 -H 2 -l)fi 2 - O 4 ] 2 + (a H 2 fl) 2 } 1/2 


(a 2 ( 


1 - n 2 ) 2 + n 2 (o 2 -i-h 2 ) 2 I' 1 


(H-7) 


where 


oi - a. I m » H t H/m 
In n 


At resonance defined by ai = ’Jo^ (0=1) , the resonant magnification 


factor is: 


RMF = [ 1+ (oH) 2 ]^ 2 /H^ 


(H-8) 


At the half-power points 


HPMF = RMF//2 


(H-9) 


Combining equations (H-7, H-8, H-9), we obtain the following 
result : 


- r-1.. 2 N 2 2 ,2 . -2. 2-,2_ r -2 ,-2 e 2 ,._2 

k la d-n hp ) + -o hp (n hp -i-H ) ] -[a -(a -H -i)n hp 


~4 -,2 , - - 2 . 2 2 

-«hp ] + < oH > % 


(H-10) 


where 


K ? [l+(i H ) 2 ] (2 H 6 )" 1 


119 



2 

Unfortunately, equation (H-10) is quartic in and thus is very 
unwieldy for practical use. 

To develop equations for the Pendered improvement of the KP method 
for a standard linear solid, we diffrentiate equation (H-6) with re- 
spect to time, obtaining the following result: 


dcp/dw = a 2 h 2 uj {[o-V - (a 2 -h 2 -uu^)uu 2 -uu 4 ] 2 + <AV } _1 (H-ll) 


At resonance defined by ufuj and denoted by subscript R, we can 
simplify equation (H-ll) as follows: 


(dcp/du)) R = |[l + (ai n /of) 2 ]uj n | 1 


or 



(H- 12) 


It is interesting to note that at resonance, the parameter h drops 
out, so that O’ can be determined directly from equation (H-12) . To 
determine h, it is necessary to measure dcp/duu at some other frequency 
(say ud') and then use equation (H-ll), with a known from equation (H-12). 
The second frequency uu' should not be too close to u;^ to avoid inaccuracy. 
However, it should not be too far from that the MF has dropped down to 
such a low value that signa 1-to-noise-ra tio difficulties develop. 

Biot Solid . - Caughey (ref. 60) presented an analysis of the sin- 
usoidally forced vibration of a single-degree-of-freedom system with 
Biot damping. In the present notation, the equation of motion is: 


120 




(H-13) 


mu + ku -g^k ^ Ei [-e(t-T)] (du/dT) dT = 1? sin &t 


where g^ and e are parameters of the model. 

The steady-state-response solution of equation (H-13) is of the 
same form as equation (H-4) , where here 


tT = (F /kO^Cl+g^tn /l+(n/u) 2 - f> 2 ] 2 +[g 1 arctan(.n/u) ] 2 j~^ 

(H— 14) 

cp = arctan (tt/ 2) g 1 [ 1+^/2) -In (1/pi) - fj 2 ] 1 (h-15) 

where 


r> = uu/ai . U 

n 


e/di 

n 


Equation (H-14) may be rewritten in terms of the magnification 
factor as follows: 


MF 


f — 2 2 

= |[l+gl'*’- n /l+(Q/u) 2 - fl ] + Lg 1 arctan (Q/u) ] 


2" 1 - 1/2 

v 

J 


( H-16) 


At resonance defined by ou = (i.e. f> = 1) , the resonant magnifi- 

cation factor is: 


RMF 


= g^ [ (-tn J 14u" 2) 2 + (arctan px" 1 ) 2 ] ^ 


(H-17) 


It is noted that the RMF defined in this way is not the peak MF, but it 
is very close to it in practical systems. Also, t is emphasized that 
the RMF is a function of both parameters, g and |j. Thus, data from some 
other frequency is required. At the half-power points (subscript hp) 
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defined by HPMF = RMF/y^ , the following relation must hold: 



Using the necessary experimental data, equations (H-17) and 
(H-18) can be used to determine the Biot material parameters g and 
U. 

To derive Pendered type KP relations, we take the derivative of 
equation (H-15) with respect to time, with the following result: 

(H - 1 9) 

At resonance defined by 0 = 1 and denoted by subscript R, we 
can simplify equation (H-19) to the following result: 

Since equation (H-20) contains both parameters (g^and y) , it is 
necessary to measure dcp/du; at some other frequency (say Q') and then 
use equations (H-19) and (H-20) to solve for the two parameters. 
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APPENDIX I 


DISCUSSION ON SYMMETRY OF THE ARRAY 
OF CONSTITUTIVE COEFFICIENTS 


In the analyses made in the body of this report, we have assumed 
that the arrays of both the elastic coefficients and the damping coef- 
ficients are symmetric, i.e. 




In the present appendix, we discuss both classical and modern think- 
ing, as well as experimental results which are pertinent to this topic 
of symmetry of the constitutive-coefficient arrays. In recent years, the 
validity of the symmetry hypothesis has been questioned for a wide 
variety of materials such as soil, geological formations, laminated wood, 
fabrics subjected to membrane loads, and certain composite materials 
(references 141-144) . 

II. Elastic Coefficients 

More than a century ago, there was considerable controversy regarding 
the symmetry of the stiffness matrix (or the compliance matrix) of a 
perfectly elastic material. About a century ago, the argument was supposedly 
resolved by the brilliant mathematical analysis due to George Green and 

123 



the experimental work of Woldemar Voigt; see reference 145* Although 
Green's analysis was purely mathematical, its importance to modern energy 
methods of static and dynamic stress analysis was recognized by Love 
(reference 146) , who stated that symmetry of the elastic coefficient array 
is a necessary condition for the existence of the strain energy function. 
Voigt verified that symmetry did hold, within experimental accuracy, for 
the wide variety of single crystals which he investigated. Thus, this 
hypothesis became well established in anisotropic elasticity theory and 
it was tacitly assumed to hold for all materials. 

Incidentally symmetry of the compliance matrix array requires that 
the following so-called reciprocal relationship must hold among four of 
the engineering elastic properties: 

= v ki^ E k < k «' t==1 * 2 * 3 ; no sum ) ( i- i) 

where the E's are Young's moduli and the v's are Poisson's ratios. 

Recently Alley and Faison (ref. 143) clearly demonstrated experimentally 
that equation (1-1) does not hold for fabrics subjected to biaxial load- 
ings. They attributed this result to fiber friction, which is a noncon- 
servative process. Table III lists some filamentary-composite-material 
test values for which eq . (1-1) does not hold. The viscosity of the epoxy 
matrix may have been the nonconservative process operative in the case 
of the glass-epoxy and boron-epoxy composites; see Section 12 of this 
Appendix. Also, the glass was in fabric form, so fabric friction could 
have played a part. Furthermore, in the case of the boron-epoxy composite, 
the minor Poisson's ratio is so very small that it is difficult to measure 
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with a high degree of confidence and experimental scatter effects become 
highly magnified. However, none of these explanations account for the 
case of graphite-graphite composite. It does not creep at the test 
temperature (room temperature); the fibers were filament wound, not 
in fabric form; and the minor Poisson's ratio is large enough to be 
measured accurately (especially since strain-gage transverse sensitivity 
was provided for in the data reduction). Thus, so far, there has not 
been a satisfactory explanation for this anomaly. 

12. Viscoelastic Material 

/ 

For viscoelastic materials which are homogeneous and anisotropic 
on a macroscopic basis, Rogers and Pipkin (ref. 141) claimed that there 
are no thermodynamic bases which require symmetry of the array of con- 
stitutive coefficients. Thus, they indicated that for such materials, 
symmetry can be established by experiment only. For composite materials 
having a viscoelastic (epoxy resin) matrix, some experimental data 
suggest that symmetry may not hold in the case of reinforcement by 
glass fibers and boron fibers (see table III). However, other experi- 
mental data on a viscoelastic composite (nylon-fiber-reinforced rubber) 
suggest that symmetry does hold. 

In view of the lack of a definite proof regarding symmetry of a 
viscoelastic material, symmetry is assumed to hold in the present in- 
vestigation (see hypothesis H6, Section 2.1). 



APPENDIX J 


COMPUTER PROGRAM DOCUMENTATION 

The program described in this appendix was programmed to ac- 
complish the following three computations: 

1. Calculate the shear factor K for laminates, using Jourawski 
static shear theory. 

2. Calculate the lowest eigenvalue for a simply supported laminated 
plate without damping. 

3. Calculate the amplitude frequency response and modified 
Kennedy-Pancu frequency response and damping data for a free- 
edge anisotropic plate with material damping. 

Computation 1 was accomplished by using an explicit algebraic expression. 
Computation 2 was performed by using IBM System/360 Scientific Subroutine 
Packages NROOT and EIGEN. Computation 3 was performed by Package SMIQ. 

A complete description of the variables, operations, etc. may be found in 
IBM Manual 360A-CM-03X, version III, for the subroutines NROOT, EIGEN and 
SMIQ. 

The program was written in FORTRAN IV language as prescribed in IBM 
System Reference Library Form C-28-6274-3. 

The input-data deck was set up as follows.: 

Computation 1 - 
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(a) Thickness of each layer (lower and upper limit) 

(b) Elastic and shear moduli 
Computation 2 - 

(a) Plate geometry 

I 

(b) Lamination geometry (specia lly orthotropic) 

(c) Moduli and Poisson's ratio data for each ply 

(d) Density for each ply 

(e) Shear factor (as calculated in Computation 1) 

Computation 3 - 

(a) Young's and shear moduli for each layer 

(b) Poisson's ratios for each layer 

(c) Bending and twisting stiffnesses for each layer 

(d) Loss tangents corresponding to moduli, Poisson's ratios, and 
stiffnesses for each layer 

(e) Plate geometry 

(f) Lamination geometry, including angle of orientation for each 
layer 

(g) Density of each layer 

(h) Shear factor (as calculated in Computation 1) 

(i) Mode numbers of the assumed modes. 

A complete listing of the computer program is presented at the end 
of this report. 
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Figure 1. Composite containing undirectiona 1 fibers. 
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Figure 2. Examples of symmetrical laminates. 
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Figure 3. An example of a balanced laminate. 





Figure 5. Lamination geometry, showing upper half of plate. '(Plane xy 
is the midplane of the plate) . 



Figure 6. Bending-stress distribution. 
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Shear factor for laminates consisting of multiple alternating layers of two materials 
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Figure 13. Lowest Eigenvalue (X) for homogeneous, isotropic 
plate. 
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Figure 19. Nodal patterns of the first five modes of composite 
material plates at angle of orientation 6 = 3C' fl 
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Figure A-l. Shingle-laminated plate. 
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SYMMETRY AXIS 


Figure A-2. Major Material-Symmetry Axis Direction. 
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THREE-PARAMETER MODELS 
(STANDARD LINEAR SOLID) 


Figure B-l. 


Two-and three -parameter viscoelastic models. 






Figure B-2 . Many-element viscoelastic models. 








(a) ELLIPTICAL HYSTERESIS LOOP 



(b) POINTED-END, STRAIGHT-SIDED HYSTERESIS LOOP 




FREQUENCY 


Figure B-4. Effect of damping on magnification factor. 
(Not drawn to scale). 



Figure B-5. Effect of damping on half-power-point frequency 
separation (U£-aj]) . (Not drawn to scale). 
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Figure B-7. Exponential decay. (Not drawn to scale). 
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N2=N2+1 

150 A(L) = A(L >*B( N1 )*X(N2 ) 



170 SUMV=SUMV+A(L)*A(L> 
175 SUMV=DSQRT(SUMV) 
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VIBRATION OF LAMINATED ANISOTROPIC PLATES 
CONSIDERING THICKNESS SHEAR AND DAMPING EFFECTS 
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i*B I26*CI 35MK*C I 55NL+B l 66* C I 36MKACI 3 7 NL ) 

JV I = JVI+ 1 

S( JJYR, JVI )=-( (A/B)*BI22*CI41MK*CI57NL+BI26*CI42MK*CI43NL+BI26 
1*C 143MKACI 4 2NLM B/A) *B I66*C I 57MK4CI 41 NL > 

JW I = J W I + 1 



S( JJYH. JWl l=-(B*K12*AI45*CI 1 5MK*CI 1 2NL+A*K22*AI 55*C 1 1 2MK4C I 1 5NL > 
JX l=JXI+i 

S< J JYR, JXI )=-{Dll2*C I25MK4C I 24NL+ I B/ A) *011 6*C 1 3 0MK*C 1 23NL+ ( A/B ) 
1*DI26*CI23MK*C130NL+DI66*CI 24MK*C I 2 5NL+ A*B*K 1 2* A1 45*C 1 2 3MK 
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INTEGRAL CI8MK.CI8NL 
IF(M-K) 1002, 1004, 1002 



1004 CI8MK=2.*M*K*PI**2 

1005 CONTINUE 

IF(N-L) 1006*1008,1006 

1006 CI8NL=0, 

GO TO 1009 

1008 CI8NL=2.*N*L*PI**2 
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INTEGRAL CI5MK.CI5NL 
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INTEGRAL CI37MK* CI37NL 
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JYR=JYR4l 

S( J JUI * JYR )=4( Ol 1 2*C I37MK*Cl 36 NL 4(B/A)*BI16 *C 15 5MK *C I 35NL4 ( A/B) 
1 *BI 26 *C I 35MK *C I 55NL4BI 66*CI 36MK*CI 3 7NL) 

JUI=JUI41 

S( J J U I » J U l )=(B/A) *AR 1 l * Cl 8MK*C I 1 NL 4 2 . * AR 1 6 *C I6M K*CI 6NL4 (A/B) 



1 * AR66 *C I IMK*CI8NL-*2*A*B*P0*CI 1 MK*C 1 1 NL 
JVI=JVI+ 1 

S(JJUUJVl »-AR12*Cl5MK.*Cl4NL+<B/A)*AfU6*CI 10MK*C13NL+1 A/Bl 
1*AR26*CI3MK*CI 1 0 NL+ AR664C I 4MK*C I 5NL 
JWI = J1*I + 1 
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IF ( N-L 1 1 406 * 14 08. 1406 
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